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' Abstract 

Integrable boundary conditions are studied for critical A-D-E and general 
graph-based lattice models of statistical mechanics. In particular, using tech- 
j> \ niques associated with the Temperley-Lieb algebra and fusion, a set of bound- 

ary Boltzmann weights which satisfies the boundary Yang-Baxter equation is 
obtained for each boundary condition. When appropriately specialized, these 
boundary weights, each of which depends on three spins, decompose into more 
natural two-spin edge weights. The specialized boundary conditions for the 
A-D-E cases are naturally in one-to-one correspondence with the confor- 
mal boundary conditions of s£(2) unitary minimal conformal field theories. 
Supported by this and further evidence, we conclude that, in the continuum 
scaling limit, the integrable boundary conditions provide realizations of the 
complete set of conformal boundary conditions in the corresponding field the- 
ories. 
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1 . Introduction and Overview 



The notion of conformal boundary conditions in conformal field theories, in the sense 
introduced in and the notion of integrable boundary conditions in integrable 
lattice models, in the sense introduced in 0, are both well developed. It is also well 
known that conformal field theories can be identified with the continuum scaling 
limit of certain critical integrable lattice models of statistical mechanics. A natural 
question to address, therefore, is whether similar associations can be made between 
the corresponding boundary conditions. We shall demonstrate here that indeed they 
can. In so doing, not only is the existence of deep connections between conformal 
and integrable boundary conditions established, but a means is also provided for 
gaining insights, generally not available within conformal field theory alone, into 
the physical nature of conformal boundary conditions. 

That there should be a relationship between such integrable and conformal 
boundary conditions is not immediately apparent and, accordingly, the correspon- 
dence is somewhat subtle. It is known that nonintegrable boundary conditions 
can be identified with certain conformal boundary conditions and, conversely, we 
expect that integrable boundary conditions will only lead to conformal boundary 
conditions upon suitable specialization. Nevertheless, we shall explicitly show for 
all of the s£(2) cases considered here that an integrable boundary condition can, 
after such specialization, be naturally associated with every conformal boundary 
condition. Furthermore, the generality of the approach used here suggests that it 
is possible to obtain an integrable boundary condition corresponding to each allow- 
able conformal boundary condition in any rational conformal field theory which is 
realizable as the continuum scaling limit of a Yang-Baxter-integrable lattice model. 

The basic context for this work is provided by certain general results on bound- 
ary conditions. For statistical mechanical lattice models, it is well known from || 
that a model is integrable on a torus by commuting transfer matrix techniques if 
its Boltzmann weights satisfy the Yang-Baxter equation. A result of || then states 
that such a model is similarly integrable on a cylinder with particular boundary 
conditions if corresponding boundary Boltzmann weights are used which satisfy a 
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boundary Yang-Baxter equation. A further result, obtained in ||, is that such 
boundary weights can be constructed using a general procedure involving the pro- 
cess of fusion. For conformal field theories, there exists a fundamental consistency 
equation of JIJ associated with the conformal boundary conditions of any rational 
theory on a cylinder. Furthermore, as shown in [|, [6], [7], || ||, the general task of 
solving this equation and completely classifying the conformal boundary conditions 
is essentially equivalent to that of classifying the representations, using matrices 
with nonnegative integer entries, of the fusion algebra. 

In this paper, we consider the critical unitary A-D-E integrable lattice models, 
as introduced in fllQ |, and the critical series of si{2) unitary minimal conformal 
field theories with central charge c < 1, as first identified in [11, and classified 



on a torus using an A-D-E scheme in |L3|, |15| . As shown in |L6|, [17], [18]], the 
continuum scaling limit of these lattice models is described by these field theories. 
Furthermore, the task of classifying the conformal boundary conditions of these 
theories on a cylinder was carried out in full in |§, leading to an A-D-E scheme 
matching that for the classification of the theories on a torus. Lattice realizations 
of some of these conformal boundary conditions have been identified and studied 
in |g, 0, gg, gg, |22], g!| [?§ |27|, ^ f2S|, but in most of these cases the 

lattice boundary conditions are not integrable in the sense of being implemented 
by boundary weights which satisfy the boundary Yang-Baxter equation. Here, we 
use the fusion procedure of H] to obtain systematically such integrable boundary 
weights and we show that, when appropriately specialized, these can be interpreted 
as providing realizations of the complete sets of conformal boundary conditions, as 
classified in 0, of the corresponding field theories on a cylinder. 

In obtaining the integrable boundary conditions for the critical unitary A-D-E 
models, it is convenient to consider a somewhat larger class of integrable lattice 



models. These models, essentially introduced in are restricted-solid-on-solid 
interaction-round-a-face models and each is based on a graph Q . In such a model, 
the spin states are the nodes of Q, there is an adjacency condition on the states 
of neighboring spins given by the edges of Q and the bulk Boltzmann weights are 
defined in terms of a particular eigenvalue and eigenvector of the adjacency matrix 
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of Q. The critical unitary A-D-E models are obtained by taking Q as an A, D or 
E Dynkin diagram and using the Perron- Frobenius eigenvalue and eigenvector of 
its adjacency matrix. These general graph-based models are also closely related to 
the Temperley-Lieb algebra, it being possible to express their bulk weights in terms 
of matrices of a certain representation of this algebra involving Q. The fact that 
these weights satisfy the Yang-Baxter equation is then a simple consequence of the 
defining relations of the algebra. 

In Section 2, we consider the abstract Temperley-Lieb algebra. The results are 
thus independent of its representation and apply to the class of lattice models as- 
sociated with the Temperley-Lieb algebra, this including the graph-based models 
of interest as well as certain vertex models. The main emphasis of Section 2 is on 
the process of fusion and its use in the construction of boundary operators which 
correspond to boundary weights. In particular, we list various important proper- 
ties satisfied by the operators which implement fusion, including projection and 
push-through properties, and we obtain several results on the properties of the con- 
structed boundary operators, including the facts that they satisfy an operator form 
of the boundary Yang-Baxter equation and that they can be expressed as a linear 
combination of fusion operators. 

In Section 3, we specialize to representations of the Temperley-Lieb algebra in- 
volving graphs and study the corresponding graph-based integrable lattice models. 
For any such model, we then obtain an integrable boundary condition and a related 
set of boundary weights for each pair (r, a) , where r is a fusion level and a is a spin 
state or node of Q. In terms of the procedure of ||], the (r, a) boundary weights are 
constructed from a fused double block of bulk weights of width r — 1, with the spin 
states on the corners of one end of the block fixed to a. These boundary weights each 
depend on three spins, but we find using results from Section 2 that, upon appro- 
priate specialization, all of the weights in a set can be simultaneously decomposed 
into physically more natural two-spin edge weights. We refer to the point at which 
this occurs as the conformal point, since it is here that we expect correspondence 
with conformal boundary conditions. While certain integrable boundary weights for 
some of the models considered here have been obtained in previous studies, specif- 
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ically ||, [7|, [H], |32], |33|, [J5|, this crucial decomposition had not previously been 
observed. We conclude Section 3 by considering in detail the symmetry properties 
of a lattice on a cylinder with particular left and right boundary conditions. 

In Section 4, we specialize to the critical unitary A-D-E models. We obtain 
completely explicit expressions for the boundary edge weights of the A and D cases 
and study A3, A4, D4 and E§ as important examples. We also obtain an important 
relation through which any A-D-E partition function can be expressed as a sum of 
certain A partition functions. For all of these A-D-E models, we find that, at the 
conformal point, the (r, a) and (g—r—1, a) boundary conditions are equivalent, where 
g is the Coxeter number of the A-D-E Dynkin diagram Q and a 1— > a is a particular 
involution of the nodes of Q. This implies that these boundary conditions are in 
one-to-one correspondence with the conformal boundary conditions of the unitary 
minimal theories A4(A g _ 2 ,Q), as classified in ||. We then use this and further 
evidence to argue that the integrable boundary conditions obtained here provide 
realizations of the Ai(A g ^2, G) conformal boundary conditions. 

In Section 5, we briefly discuss ways in which the formalism of this paper could 
be applied to other models. 



2 . Relevant Results on the Temperley-Lieb 
Algebra 

In this section, we list and obtain various results on the Temperley-Lieb algebra. 
The defining relations of this algebra were first identified in |36]| and the formalism 
used here is largely based on that of |37|, [39|| . 

Our primary objective in this section is to study operators in the abstract algebra, 
which, in certain representations of the algebra, correspond to bulk weights and 
boundary weights of a lattice model. In particular, motivated by the procedure of 
||, we shall construct certain boundary operators. This procedure involves fusion, 
which can be regarded as a means whereby new fused operators satisfying the Yang- 
Baxter equation are formed by applying certain fusion operators to blocks of unfused 
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operators. The process of fusion was introduced in ||0 and first used in the context 
of the Temperley-Lieb algebra in |37|, |38 . 



2.1 The Temperley-Lieb Algebra and General Notation 

The Temperley-Lieb algebra T(L, A), with L G Z>o and A G C, is generated by the 
identity / together with operators e±, . . . , ex which satisfy 

= 2 cosA tj 

tj e k ej = ej, \j-k\ = 1 (2.1) 
ej e k = e k e j} \j-k\ > 1 . 

The various operators to be studied in Section 2 will all be elements of T(L, A) 
for some fixed L and A. 

Throughout this paper, we shall use the notation 

f Sin(rA + M) , A/vr^Z 
s r ( u ) = I smA (2.2) 

[(_l)(r+l)A/*r (r +u ) } A/ttGZ, 

for any r G Z and u G C, with A being the same constant as in T(L, A). 

When T(L, A) is applied to a lattice model, A is the crossing parameter. For the 
A-D-E models, this parameter is always a rational but noninteger multiple of 7r 
so that the first case of ( |2.2| ) applies. We also note that the second case of ( |2.2|) is 
simply a limiting case of the first, 

, , , sinfrA'+fA'-AW) , , 

sJu) = lim — i L LI 2.3 

v ' A/^eZ a'-a sinA' ' V ' 

so that in proving identities satisfied by these functions it is often sufficient to 
consider only the first case. 

We shall also denote, for any r G Z, 

S r = s r (0) . (2.4) 

We see, as examples, that for any A, 

S = 0, Si = 1, 5 2 = 2cosA. (2.5) 
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2.2 Face Operators 

We now introduce face operators Xj(u), for j E {1, . . . , L} and u G C, as 

Xj(u) = si(—u) I + Sq(u) e 3 . (2.6) 

These operators correspond in a lattice model to the bulk Boltzmann weights of 
faces of the lattice and in this context u is the spectral parameter. 
We shall represent the face operators diagrammatically as 



X 3 {u) 




(2.7) 



J-l 3 3+1 

From the Temperley-Lieb relations ( |2.1| ) and properties of the functions ( |2.2| ), it 
follows that the face operators satisfy the operator form of the Yang-Baxter equa- 
tion, 

X 3 (u) X j+1 (u+v) X 3 {v) = X j+1 (v) Xjiu+v) X j+1 (u) 



(2.8) 



\ v 














( u+v \{ 


\ u 








3-1 3 


3+1 3 


+2 j- 


-1 3 3+1 3+2 



We see that the face operators also satisfy the commutation relation 



X 3 iv)X 3 iu) 



and the operator form of the inversion relation 

Xj(—u) Xj(u) = Si(—u) Si(u) I . 



(2.9) 



(2.10) 



2.3 Fusion Operators 

We now proceed to a consideration of some aspects of the process of fusion. In 
particular, we shall state some important properties of fusion operators. The proofs 
of these or similar properties can be found in |37], . 
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In general, fusion levels correspond to representations of a Lie algebra. In the 
case considered here, this algebra is s£(2) and the fusion levels are labeled by a single 
integer r & {1,2, ... , g}, where the maximum fusion level g depends on A according 
to 

h, A = kn/h, k and h coprime integers, h > 1 
oo, otherwise . 
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(2.11) 



We note that for the A-D-E lattice models, the first case of (|2.11|) always applies. 

We now introduce fusion operators PJ, for r G {1, . . . , min(g, L + 2)} and j G 
{!,..., L + 3 — r}, these being defined recursively by 



P 1 = P 2 = I 

3 3 

P] = ^- P£l X,(-(r-2)A) P^l , r > 3 . 



f2.12l 



We note that, for finite g, the restriction of fusion levels to r < g is necessary in 
order to avoid SV-i — in the denominator in (|2.12|) . 

We shall represent the fusion operators diagrammatically as 



(2.13) 



j-l j j+r-3 j+r-2 

In general, PJ can be expressed in terms of / and ej, . . . , ej +r _3, the cases of the 
first few fusion levels being 



P 1 = P 2 

3 3 



1 



1 ' ^ 



(2.14) 



Pj = I -—(e j + e j+1 ) + —(e j e j+1 + e j+1 e j ) 



S; 



s* 
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The fusion operators can also be expressed as 




j-l j j+r-3 j+r-2 j-1 j j+r-3 j+r-2 



A key property of the fusion operators is that they are projectors, 

(pry = pr (2J6) 

In fact, more generally, we have the property 

Pj! PJ = PJ Pj! = PJ , if < f-j <r-r'. (2.17) 
A useful case of this is r' = 3, from which we have 

e r PJ = PJe r = 0, ifj<j'< j+r-3. (2.18) 
The fusion operators also satisfy the Hecke-like identity 

C<2 pr pr pr pr _ q2 pr pr pr pr _ q q pr+1 icy i q\ 

2.4 Fused Row Operators 

We now introduce operators YJ(u) and YJ"(u), for r G {1, . . . , min(g — 1, L+l)} and 
j G {1, . . . , L+2 — r}, which correspond to fused rows of faces. These operators are 
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related to products of r — 1 face operators, for r > 2, by 

= Pr +1 X 3 (u-(r-2)\) . . . X j+r ^{u-X) X j+r ^(u) 
- Xj(u)X j+1 (u-X)...X j+r _ 2 (u-(r-2)X)PT 



r f[s k (-u) YJ 

k=i 



u 





and 



Us-k(-u) YJ(u) 

k=0 



j-1 j j+r-2 j+r—1 j—\ j j+r-2 j+r—1 

Pj X j+r _ 2 (u) X j+r _ 3 (u+X) . . . X>+(r-2)A) 
X j+r _ 2 (u + (r-2)X) . . . X j+1 (u+X) Xj(u) PJ +1 





(2.20a) 



(2.20b) 



3-1 j j+r-2 j+r-l j-1 j j+r-2 j+r-1 

The second (or fourth) equalities in (2.20) imply the push-through relations 

YJ{u) = P; +1 YJ(u) = YJ-(u)PJ 

YJ(u) = P^YJin) = YT(u)PT +1 , 
and can be derived using ( |2.15| ) and repeated application of (pT8|) . 

The fused row operators can be written in terms of fusion operators as 

YJ(U) = S r „ lSl ( U )PJ +1 PJ - SrSo^PJ^ 

Yf(u) = S r Aa)PJ - S r s r _ 2 (u) PJ +1 . 



(2.21) 



(2.22) 
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The equivalence of (2.20) and ( |2.22| ) follows from several of the properties of fusion 
operators listed in Section f2~3|. 



We note, as examples, that for the first two fusion levels, 

Y}{u) = -s {u)I, Y}{u) = s x {-u)I, Y?(u) = Y?(u) = Xj(u) . (2.23) 

An important property of the fused row operators is that they satisfy the mixed 
Yang-Baxter equations, 

Xj(u-v) Yj +1 (u) YJ(v) = YJ +1 (v) YT(u) X j+r ^(u-v) (2.24) 

X j+r _ x {u-v) YJ{u) YJ +1 (v) = YJ(v)YJ +1 (u)X,(u-v) (2.25) 

YJ(u) X j+r _ x {u+v) YJ(v) = YJ +1 (v)X J (u+v)YJ +1 (u). (2.26) 



These equations can be obtained using (2.20) and repeated application of (|2.8| ). 
The fused row operators also satisfy the product identities 

YT( u )YJ{v) = Sl (u) s r ^(v) PJ +1 - S r s (u+v) PJ +1 

(2.27) 

Y((u)Yf(v) = s r ^(u) Sl (v) PJ - S r s (u+v) P[ +1 , 



these being most easily derived using ( |2.22|) and properties of the fusion operators. 



2.5 Boundary Operators 

We now introduce boundary operators i^J(u,^), with ( G C, as products of two 
fused row operators, 

Kj(u,C) = -Yf(u-\-0YJ(u + \+0- (2.28) 

These operators correspond in a lattice model to boundary Boltzmann weights and 
in this context £ is a boundary field parameter. 

From ( p.27| ), we see that the boundary operators can be written in terms of 



fusion operators as 

iq(u,0 = s (£-u) s r (Z+u) PJ +1 + S r s (2u) PJ +1 . (2.29) 

11 



An alternative expression, which follows using ( |2.12| ) on the second term on the right 
side of ( ggg ), is 

kj(«>0 = s (z+u)s r (t-u)pj +l + ^=i So (e-«) Sr .(e+ w )pj; ie ,p; +1 . (2.30) 

We see, as examples, that for the first two fusion levels, 

Kj(u,0 = so(Z+u) Sl (e-u) I, K](u,0 = so(Z+u) s 2 (£-u) / - s (2 M ) e,, (2.31) 

The key property of the face and boundary operators is that they satisfy the 
operator form of the boundary Yang-Baxter equation, 



X^u-v) Kr +1 (u,0 Xj(u+v) Kr +1 (v,£) = 

K] +1 (v,£) Xfa+v) K] +1 (u,0 Xjiu-v) 



(2.32) 



This can be proved by substituting (|2.28|) into the left side of fl2.32|) , using (|2.26| ) 
followed by ( |2.24|) to bring Xj +r _i(u—v) adjacent to Xj +r _i(u+v ), interchanging the 
order of these face operators using (|2.9| ), and then using (|2.25|) followed by (|2.26|) to 



give the right side of (|2.32|) . 

In terms of the construction procedure of Q, the boundary operators Kj(u, £) 
can be considered as a family of solutions, one solution for each value of r, of (|2.32|) 
with given Xj(u), these solutions having been constructed by starting with the 
identity solution of ( p.32[) and adding a fused double row of faces of width r — 1. 

The boundary operators also satisfy the operator form of the boundary inversion 
relation, 

#J(u,0 = s (Z-u)s (£+u)s r (£-u)s r (Z+u)P; +1 , (2.33) 

this being most easily obtained using ( ^.29[ ) and properties of the fusion operators. 
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3 . Lattice Models Based on Graphs 



In this section, we consider representations of the Temperley-Lieb algebra involving 
a graph and the associated graph-based lattice models. Our general treatment is 
motivated by that introduced in (TO, BTJ and, when considering fusion for these 



models, our approach is also motivated by certain results of |^TJ, ^2], ^3], [44 1. 

A graph-based lattice model of the type to be considered here can be associated 
with any pair Q and ip, where Q is a connected graph containing only unoriented, 
single edges and ip is an eigenvector of the adjacency matrix of Q with all nonzero 
entries. All of the formalism and results of Section 3 are applicable to any choice 
of Q and tp, and only in Section 4 will we eventually specialize Q to be an A, D or 
E Dynkin diagram and ip to be the Perron- Frobenius eigenvector of its adjacency 
matrix, thus giving the critical unitary A-D-E models. 

In a lattice model based on Q, a spin is attached to each site of a two-dimensional 
square lattice, with the possible states of each spin being the nodes of Q and there 
being a lattice adjacency condition stipulating that, in any assignment of spin states 
to the lattice, the states on each pair of nearest-neighbor sites must correspond to 
an edge of Q. These models are interaction-round-a-face models, so that a bulk 
Boltzmann weight is associated with each set of four spin states adjacent around a 
square face. Here, we shall also obtain and use sets of boundary weights, each of these 
weights being associated with three adjacent spin states. The partition function of 
the model on a cylinder is then the sum, over all possible spin assignments, of 
products of Boltzmann weights, with each square face in the bulk of the lattice 
contributing a bulk weight and each alternate triplet of neighboring sites on the 
boundaries contributing a boundary weight. 

The key property of the boundary weights obtained in this section is that they 
satisfy the boundary Yang-Baxter equation for interaction-round-a-face models. 



This equation was first used in [33, 45] and is based on the reflection equation intro- 
duced in PEfl . However, we note that the boundary weights obtained here depend on 
certain additional fusion indices, which gives them and the boundary Yang-Baxter 
equation they satisfy a somewhat more general form than the form of those used in 
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all previous studies. We shall show that a further important property of the three- 
spin boundary weights used here is that, at a certain point, they simultaneously 
decompose into more natural two-spin boundary edge weights. 

Finally, using these and other local properties of the bulk and boundary weights, 
we shall identify various symmetry properties of the partition function and transfer 
matrices, including the invariance of the partition function under interchanging parts 
or all of the left and right boundaries, and the facts that the transfer matrices form 
a commuting family and, in certain cases, have all nonnegative eigenvalues. 

3.1 Graphs and Paths 

Throughout Section 3, we shall be considering a finite graph Q with an associated 
adjacency matrix G. We require that Q contain only unoriented, single edges, im- 
plying that G is symmetric and that each of its nonzero entries is 1. In fact, the 
formalism of Section 3 can be generalized straightforwardly to encompass graphs 
with multiple edges, but since this is not needed for the A-D-E cases of primary 
interest, we shall for simplicity restrict our attention to graphs with only single 
edges. 

We also require that Q be connected, implying that G is indecomposable. Perron- 
Frobenius theory then implies that G has a unique maximum eigenvalue with an 
associated eigenvector whose entries are all positive. 

We shall denote the set of all r-point paths on Q, for r e Z>i, by Q r ; that is, 



We note that Q l corresponds to the set of nodes of Q and that Q 2 is the set of edges 



We shall also denote the set of all r-point paths between a and b in Q by Q r ah \ 
that is, 



r-2 




(3.1) 



of g. 




(3.2) 



It follows that 



(G r -%. 



(3.3) 
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3.2 Graph Representations of the Temperley-Lieb Algebra 

A graph representation involving Q of the Temperley-Lieb algebra T(L, A) exists for 
each A for which 2 cosA is an eigenvalue of G with an associated eigenvector if) whose 
entries are all nonzero; that is, for each A for which there exists a vector if> satisfying 

G a bipb = 2cosA ipa and ip a ^ , for each a G Q. (3.4) 

be Q 

There is always at least one such case, namely that provided by the Perron-Frobenius 
eigenvector and eigenvalue. We shall assume, throughout the rest of Section 3, that 
a fixed choice of A and ip has been made. We note that, since G is symmetric, cosA 
must be real so that each if> a can, and will, be assumed to be real. In the rest of this 
paper, we shall often use ipf and ip^ 4 , which we take to be positive for positive if) a 
or to have arguments 7r/2 and 7r/4 respectively for negative if> a . 

The elements of the representation of T(L, A) associated with Q and if) are ma- 
trices with rows and columns labeled by the paths of Q L+2 , with the generators ej 
being defined by 

<^_ia i+ i II 5a k b k ■ (3-5) 



e 



3 (ao,...,a L+1 ),(b ,...,b L+1 ) ' , 



It follows straightforwardly that the defining relations ( |2.1|) of T(L, A) are satisfied, 
with the first relation depending on 

3.3 Bulk Weights 

We now proceed to a consideration of the lattice model based on the graph Q and 
associated with the adjacency matrix eigenvector if}. 

The bulk weights for this model are given explicitly, for each (a, b, c, d, a) e Q 5 , 

by 

W ( d \ < \ & S -M ^ fL r 
la b J = Sl (~ u ) ° ac ^ ~1 " M ' 

where u is the spectral parameter and A, on which s depends through ( |2.2| ), is 
the crossing parameter. The spectral parameter can be considered as a measure 
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of anisotropy, with u = A/2 being an isotropic point and u = and u = A being 
completely anisotropic points. We note that the number of bulk weights is tr(G 4 ). 
We shall represent the bulk weights diagrammatically as 



W 



d 


c 














•) 






u 


[a 


b 




a 


\ 



(3.7) 



These bulk weights are related to the face operators 
d c 



by 



W 



a b 



u 



^'( w )(eo,...,ej_2,d,o,6,ej +2 ,...,e£ + i),(eo,...,ej_2,rf,c,6,ej + 2,...,ei + i) > (3- 



where Xj(u) is taken in the graph representation of T(L, A). 
We see that the bulk weights satisfy reflection symmetry, 



d c 
a b 



u 



W\ 



d a 
c b 



it 



b c 
a d 



u 



crossing symmetry, 



W 



and the anisotropy property 





c 








•) 




b 







d 




(: 






a 





W\ 



d c 
a b 







Sac ■ 



(3.9) 



(3.10) 



(3.11) 



It follows from (|2.8|) , ( |2.10 ) and (|3.8| ) that the bulk weights also satisfy the Yang- 
Baxter equation, 

If e 



see 

( G bg G dg G fg = V 
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"{I 


d 


•) 


\ a 


b 


u — vj 




c 





aeG \ 



g&5 

{GagGcg Geg — 1) 



/. . ./ 



a g 
b c 



v W\ 



f e 
a g 



u) W 



e d 
9 c 



u — v 



(3.12) 




V 

\ 


d = a 


u 

\ , 






u 




V 


\ 




\ 




b b 



16 



for each (a, b, c, d, e, /, a) G Q 7 , and the inversion relation 



d d 




(3.13) 



= siiu) si(-u) <5 ac , 



for each (a, 6, c, d, a) G (? 5 . In these and all subsequent diagrams, solid circles are 
used to indicate spins whose states are summed over and dotted lines are used to 
connect identical spins. 



3.4 Fusion 

We now consider various objects related to the process of fusion in graph-based 
lattice models. 



3.4.1 Maximum Fusion Level 

Having chosen a graph Q and an adjacency matrix eigenvector with eigenvalue 2 cosA, 
the maximum fusion level g is then determined by A according to (|2.11| ). We note 



that although A itself is only determined by the eigenvector up to sign and shifts of 
2tt, such changes in A do not affect g or any other properties of interest. 

As we shall see in more detail in Section 4.1, if Q is an A, D or E Dynkin 
diagram, then g is the Coxeter number of Q, since any eigenvector of G with all 
nonzero entries has an eigenvalue 2 cos(kir/h), where h G Z>2 is the Coxeter number 
of Q and k G {1, . . . , h— 1} is a Coxeter exponent coprime to h. 

Although the A-D-E cases are of primary interest, simply as a further example, 
we briefly discuss the A^\ D^ l > and E^ Dynkin diagrams of the affine Lie algebras. 
If Q is one of these graphs, then any eigenvalue of G can be written as 2 cos(k7r/h), 
where h is the Coxeter number of Q and k G {0, . . . , h} is a Coxeter exponent. For 
all of these graphs, k = is a Coxeter exponent, giving a maximal eigenvalue of 2, 
and in some cases k = h is also a Coxeter exponent, giving a minimal eigenvalue of 
—2, so if either of these eigenvalues is chosen, then g = oo (and, incidentally, the 
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second case of ( |2.2| ) applies). However, if an eigenvalue corresponding to a Coxeter 
exponent < k < h is chosen, then g is finite with g < h. 

3.4.2 Fusion Matrices 

We now introduce, for each r e {1, . . .,g} and a,b G Q satisfying {G r ~ l ) a b > 0, a 
fusion matrix P r (a, b) with rows and columns labeled by the paths of Q r ab and entries 
given by 

P T (a, &)(a,ci,...,£v- 2 ,&),(a,di,..,d r ^2,&) = 

'1, r = l (3.14) 

-^T( a ,ci,...,c r _2,fe),(a,d lv ..,d r _2,6) ' T = 2 ' " " " ' 9 i 



where the fusion operator P{ is defined by (|2.12|) and taken in the graph represen- 
tation of T(r — 2, A). 



It follows from ( |2.16| ) that each fusion matrix is a projector, 

P r {a,b) 2 = P r {a,b). (3.15) 



It also follows, from fl2.12| ), ( |3.8| ) and the first equality of (|3.9|), that each fusion 
matrix is symmetric, 

P r (a,b) T = P r (a,b) (3.16) 



and, from (|2.15|) and (p.8|) , that a and b can interchanged according to 



P r (a, &)(a,ci,..., Cr _2,fc),(a,di,...,d r _2,fe) - P r (b, Ct)(fe, Cr _ 2 ,..., Cl ,a),(Mr-2,...,rfi,a) • ( 3 -17) 

We note that the fusion matrices also satisfy more general projection- type relations 
corresponding to those, (|2.17|) , satisfied by the fusion operators. 

3.4.3 Fused Adjacency Matrices 

We next introduce fused adjacency matrices, F 1 , . . . , F 9 , which are defined by the 
s£(2) recursion, 

F 1 = J; F 2 = G; F r = G F r ~ 1 — F r ~ 2 , r = 3, ...,g. (3.18) 
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We see that each fused adjacency matrix is a polynomial in G and hence that the 
set of these matrices is mutually commuting. In fact, the polynomial form is given 
by the Type II Chebyshev polynomials U r , 

F r = U r -x{G/2). (3.19) 

We note that at various places in Section 4, we shall be considering the fused adja- 
cency matrices of two different graphs and in these cases we shall explicitly indicate 
the dependence on the graph as F{Q) r . 

The main relevance of the fused adjacency matrices at this point is that their 
entries give the ranks of the fusion matrices, according to 

fO, (G p - 1 U = 
Kb = (3-20) 
(rankP'-M), (G r ~ 1 ) ab > . 

This result can be proved by defining matrices F r with entries F£ b given by the 
right side of ( |3.20|) and showing that these matrices satisfy relations (|3.18|) , so that 



F r = F r . Showing that F r satisfy the initial conditions of (|3.18 ) is straightforward. 



Meanwhile, showing that F r satisfy the recursion relation of ( 3.18Q can be done by 



using the fact that the rank of any projector is given by its trace (since by idempo- 
tence each eigenvalue is either or 1, the rank is the number of eigenvalues which 
are 1 and the trace is the sum of eigenvalues) and then using definitions ( |2.6|) , ( |3.5|) 
and Q3.14 ), the recursion relation of (|2.12| ), which is needed twice, relations fl2.16p 



and (|3.4|), and general properties of the trace and of the numbers S r . 

We note that ( |3.2U| ) also implies that F r ah are nonnegative integers. In fact, the 
property that the entries of F r are integers follows immediately from the properties 
that the entries of G are integers and that each F r is a polynomial, with integer 
coefficients, in G, but the property that the entries of F r are also nonnegative is less 
trivial and depends on the restriction r < g and on some of the additionally-assumed 
properties of G, such as its being symmetric. Since the result of nonnegativity is 
of some interest in its own right, we note that ( |3.20| ) specifically implies that if 
G is any symmetric indecomposable matrix with each entry in {0, 1} and if F r is 
defined in terms of G by ( |3.19|) , then each entry of F r is a nonnegative integer for 
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any r G {1, . . . , g m ax}, where g max is the largest of the values g corresponding to the 
eigenvectors of G with all nonzero entries. 

3.4.4 Fusion Vectors 

It can be shown using standard results of linear algebra that any symmetric idem- 
potent matrix P = P T = P 2 with complex entries can be orthonormally diagonal- 
ized; that is, written as P = UDU^ 1 , where the matrix of eigenvectors U satisfies 
U^ 1 = U T and D has rankP l's on the diagonal and all other entries 0. We note 
that if P = P* then this simply amounts to an orthonormal diagonalization of a real 
symmetric matrix, but that if P ^ P* then it differs from a unitary diagonalization 
and that, in fact, P is then not normal and a unitary diagonalization is not possible. 

From this result, it follows that each fusion matrix P r (a, b) can be decomposed 
using F£ b orthonormal eigenvectors with eigenvalue 1. We shall denote, for r G 
{l,...,g} and a,b G Q satisfying F r ah > 0, such eigenvectors as U r (a,b) a , where 
a — 1, • • • , F^ b , and refer to these as fusion vectors. 

The decomposition and orthonormality are then 

F ab 

J2 U r (a,b) a U r (a,b) a T = P r (a,b) (3.21) 

a=l 

and 

U r (a,b) a T U r (a,b) a , = <W, (3.22) 

or, more explicitly, 

ao 

^ U { a , &)a,(a,ci,...,c T w 2 ,6) U (a, fc) a ,(a,di ,...,d r -2,b) = P ( a > ^) (a,ci,...,c^-2,b),(a,di,...,d.^J2,b) 
a=l 

and 

U r (a,b) a ,(a,c u ...,c r - 2 ,b) U r (a,b) a >^ Cu ..^ Cr _^ b ) = 8 aa > . 

(a,a,...,c r -2,b) G Q r ab 

We note, as examples, that for the first two fusion levels, 

P 1 (a,a) {a)M = ±U 1 (a,a) h ( a ) = 1, for each a G Q , (3.23) 
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and 



P 2 (a, 6)( ,6),(a,6) = ±^ 2 (a,&)i,(a,fe) = 1, for each (a,b) G Q 2 . (3.24) 

We shall assume that a specific choice of U r (a,b) a has been made. All other 
possible choices are then given by transformations, 

Kb 

U r (a,b) a i-> J2 R r (a,b) aa ,U r (a,b) a ,, (3.25) 

a'=l 

where R r (a,b) is an orthonormal matrix, 

ab 

R r (a,b) aa/I R r (a,b) a , a// = 5 aa > . (3.26) 

a"=l 

We shall see that all of the lattice model properties of interest will be invariant under 
such transformations and thus independent of the choice of fusion vectors. 



3.5 Boundary Weights 

We now use the boundary operators Q2.28|) to obtain a set of boundary Boltzmann 
weights for the lattice model for each pair (r, a), where r G {1, . . . , g — l}isa fusion 
level and a is a node of Q. It is thus natural to regard these pairs as labeling the 
boundary conditions, 

{boundary conditionsj < — ► {(r, a) \ r G {1, . . . ,g— 1}, a G Q\ . (3.27) 

The (r, a) boundary weights are given, for each (6, c, d) G Q 3 with F£ a F^ a > and 

pe{i,...,F b r a },5e{l,...,F% a },by 



TT^rm 1 1/4 U r (b,a)( 3 ^ eu ..^ r _ 2)a) U r (d } a) S) (d,f 1 ,...,f r ^,a) 

SoWPb ( fc , ei ,..., er _ 2 , a)e g L 
(dj u ...j r - 2 ,a)eg r da 



,e r -2,a),(c,d,fi,---Jr-2,a) J 
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where K[(u, £) is taken in the graph representation of T(r— 1, A) and £ corresponds 
to a boundary field. We shall represent the boundary weights diagrammatically as 

d 8 



B ra c 



d 5 

b (3 



«,£ 




(3.29) 



b (3 



We see that 



number of (r, a) boundary weights = ((F r ) G 



(3.30) 



We shall refer to a boundary weight (|3.28 ) as being of diagonal type if b — d 
and (3 = 5 and as being of nondiagonal type otherwise, and we note that in the 
majority of previous studies involving such boundary weights, only diagonal weights 
were considered. 

We also note, as an example, that the (1, a) boundary weights are all diagonal 
and given by 



B la [c 



a 1 
a 1 



(3.3i; 



so(2£) W 

A (1, a) boundary condition is thus one in which the state of every alternate bound- 
ary spin is fixed to be a, while each other boundary spin, whose state can be any c 
adjacent to a on Q, is associated with a weight proportional to ipf . We shall refer 
to such a boundary condition, in which the state of every alternate boundary spin 
is fixed, as semi-fixed. 

We now find that the boundary weights can also be expressed as 

r- 2 



d s 

b (3 



u,A = - (s (2£) J] s-fc-i(«-0 Sk(u+0) 1 



k=l 



U r (d } a) St{d j 1 j 2t ...j r _ 2t a) 

fl h Sr-2 





-(r-2)A 

\ 


-«-{ 
-(r-3)A 

\ 




-u-i 

\ 


WW " 




u-i 








-(r-l)A 


-(r-2)A 






\ 


\ 




\ 



(3.32) 



U r (b,a 



/3,(fe,ei,e 2 ,...,e r _2,a) 
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this following straightforwardly from ( ggg ), (2.20), ( gTg ), ([T6|), ( |3~T0D , flg^IQ and 
( |3.22| ). In this form, we see that the (r, a) boundary weights can be considered as 
having been constructed by starting on the right at node a, essentially with the level 
1 boundary weights ( |3.31| ), and adding a fused double row of bulk weights of width 
r — 1. 



The boundary weights, together with the bulk weights (|3T6| ), satisfy the boundary 
Yang-Baxter equation, 

(g,h) e G 2 7=1 V 
(G bg G dg G eh F^ a >0) 




na 



E E Br 

(g,h)eg 2 7=1 
(G ch G dg G /s F- a >0) 



ft, 7 
6 /3 



v,aw< 



g h 
d c 



X—u—v x 



B ra \g 



f 5 
h 7 



f 9 

e d 



(3.33) 



u — v 




f / 




for each (6, c, d, e, /) G £ 5 with FJ a > and /3 G {1, . . . , F b r a }, 5 G {1, . . . , FJJ. 
This equation can be verified by expressing the bulk and boundary weights in terms 
of bulk and boundary operators, using (|3.8|) and (|3.28|) , and applying the operator 



form of the boundary Yang-Baxter equation, (|2.32|) . In doing this, the fusion matrix 
which is formed on the interior of each side of the equation by the sum on 7 and 
Q3.21 ) can be moved to the exterior by expressing the fusion matrices in terms of 
fusion operators using ( |3.14 ) and the boundary operators in terms of fused row 
operators using ( |2.2£j| ) and applying the push-through relations ( |2.21| ). Meanwhile, 
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the orientation of the central bulk weights on each side of ( |3.33| ) can be changed 
to that required for ( p.32j ) by using crossing symmetry fl3.1Q|) , with the eigenvector 
entries from ( |3.10| ) canceling where needed with those introduced explicitly through 

(P3BD - 

We note that the boundary Yang-Baxter equation (|3.33|) is still satisfied after 
renormalization of the boundary weights, 



B ra c 



d 5 

b (3 



uA -» f ra (u,OB ra (c d b 



d 5 

P 



(3.34) 



where f ra are arbitrary functions. It is also still satisfied after gauge transformations 
of the boundary weights, 



B ra \ c 



d 5 

b (3 



«,( ^ E E S ra {b) m S ra (d) ss , B™ [c 



d 6' 

b f3' 



where 5" ra (e) are arbitrary orthonormal matrices, 



a"=l 



S ra (e) aa i> S ra (e) a / a /i — 5 aa i 



(3.35) 



(3.36) 



Indeed, a transformation ( p.25|) of the fusion vectors simply induces a gauge trans- 
formation of the boundary weights with 



S ra ( e ) = R r (e,a) 



(3.37) 



In addition to the boundary Yang-Baxter equation, some other important local 
relations satisfied by the boundary weights are boundary reflection symmetry, 



B ra \ c 



d 5 

b (3 



B ra \ c 



b (3 
d 5 



(3.38) 



boundary crossing symmetry, 



E w( 



(b,e,d) e G 



d e 
c b 



2m- A B ra e 



d 6 

b (3 



s (2u) B ra \c 



d 5 

b (3 



(3.39) 
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and the boundary inversion relation, 

F r 



ea I V4 1 .1/4 j 1/2 , 

yfe re £ro/ 



e£5 7=1 

(G ce F e r a >0) 



e 7 
b (3 



u, £ B ra c 



d 5 

e 7 



s (£—u) s (£+«) s r (C-w) s r (6+«) 



(3.40) 



8bd S/3S 



for each (b, c, d) G G 3 with i^i^ > and (3 G {1, . . . , F b r a }, 5 G {1,...,FJJ. 
Boundary reflection symmetry follows from (|3.28|) , boundary crossing symmetry 
can be proved by expressing the boundary weights in ( |3.39|) in the form (|3.32|) and 
repeatedly applying the Yang-Baxter equation ( |3.12| ), and the boundary inversion 
relation follows from the operator form of the boundary inversion relation, 



3.6 Boundary Edge Weights 

We now introduce boundary edge weights in terms of which the boundary weights 
of the previous section can be expressed. 

We begin by defining, for each boundary condition (r, a) , a set of boundary edges 

as 

£ ra = {(b,c)eg 2 \F£ a F r ca +1 >0}. (3.41) 

We note that the boundary edges are ordered pairs and that, in contrast to the 
graph's set of edges for which (6, c) G Q 2 <^ (c, b) G Q 2 , the appearance of (6, c) in 
£ ra need not imply the appearance of (c, b) in S ra . 

The (r, a) boundary edge weights are now given, for each (b, c) G S ra , (3 G 
{1,. ..,!£} and 7 e{l,...,F;+ 1 },by 

E ra (b,c) Pl = 



S?W sr nrf h , rrr+ lr \ (3-42) 

(b,di,...,d r -2,a) S 5J a 



V'6 



We shall represent the boundary edge weights diagrammatically as 



E ra (b,c)^ = . (3-43) 

b p 
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We see that 



number of (r, a) boundary edge weights = (F r G F 



(3.44) 



We note, as an example, that for the (l,a) boundary condition, 
S la = {(a,c) | G ac = 1}, E la (a,c) n = ±^f/^f. 



(3.45) 



We now find, substituting (^2|) into ( ggg ) and using ( gig ), ( ggg ) and ( ggg) , 
a general expression for the boundary weights in terms of boundary edge weights, 



B ra c 



d 5 



6 /3 



so(20 <A 



1/2 



F r+1 

*w*/» + £ 4S £ E ra (b,c) M E ra (d,c) 5l 



(3.46) 



From this, we immediately see that at u = £ the boundary weights are indepen- 
dent of £ and can be decomposed as a sum of products of boundary edge weights, 



B ra c 



d 5 

b (3 



/ 7=1 

r, a 





(3.47) 



r, a 



We note that the origin of this decomposition is the fact, apparent from 02.29 ), 
that at it = ^ the boundary operators are proportional to fusion operators, so that 
the decomposition of boundary weights is essentially equivalent to the eigenvector 
decomposition of projectors. 

We shall refer to this point, u = £, as the conformal point, since it is here that 
certain lattice models are expected to exhibit conformal behavior, with the set of 
(r, a) boundary edge weights providing a lattice realization of a particular conformal 
boundary condition. 
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We find, by using ( |2.30| ) in ( |3.28| ), that a decomposition similar to Q3.47D also 



occurs at u 



B ra \ c 



d 5 

b (3 



75 



(3.48) 



7=1 

(rjfcl) 



We therefore see that, apart from an unimportant reversal of the order of the nodes 
in the boundary edges, the (r, a) boundary condition at u = — £ is equivalent to the 
(r — 1, a) boundary condition at u = £. In fact, decompositions of the form ( |3.47| ) 
or ( |3.48[ ) also occur at other points, for example at u = — rA— £, u = rA+£ or points 
related to these by trigonometric periodicity, but since these points all involve the 
same boundary edge weights up to reordering of the nodes in the boundary edges 
or relabeling of r as r±l, we can, without loss of generality, restrict our attention 
to the single conformal point u = £. 

We also note that the sum in (|3.47| ) is empty if -FJ a +1 = 0. It is thus possible that 
for a particular boundary condition, certain boundary weights are non-zero away 
from the conformal point but vanish at the conformal point. From (|3.46|) , we find 
that such boundary weights are specifically those for which b = d, j3 = 5, F£ a Gb c > 
and F? 1 = 0. 

We see from (|3.46|) that the boundary weights satisfy the boundary anisotropy 
property 



B' 



d 5 

b (3 



s (0 *r(0 V? 



s (2£) K 

The second term on the right side of ( |3.46| ) also vanishes for £ 
assuming X/n (jL Z, 



d 5 

b (3 



it, ±iooj 



± 



1 e 



2 sinA ipb* 



8bd 5/36 



(3.49) 
±ioo so that, 

(3.50) 



Comparing the right sides of ( |3.49 ) or (|3.50|) with the (1, a) boundary weights ( p.31|) , 
we see that at u = or £ — > ±ioo the nonzero (r, a) boundary weights reduce, up 
to unimportant normalization, to (1,6) boundary weights, with each (1,6) weight, 
for any 6, appearing F£ b (which may be zero) times. 
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Finally, we note that a gauge transformation of the boundary edge weights, 

FT F r+1 

oa 1 ca 

E ra (b,c)^ -> £ E S ra (b)wS r+1 ' a (c) 77/ E ra (b,c) 0lj ,, (3.51) 

/3'=1 7 '=1 



induces a gauge transformation ( [3.35] ) of the boundary weights, for any S ra (e) sat- 



isfying (|3.36|) 



3.7 Transfer Matrices and the Partition Function 

We now proceed to a study of some aspects of the complete lattice. We shall be 
considering a square lattice on a cylinder of width N and circumference 2M lattice 
spacings, with the left boundary condition and boundary field given by (r 1; ai) and 
£1 and the right boundary condition and boundary field given by (r 2 , a 2 ) and £ 2 . In 
particular, we shall express the partition function for the model using double-row 
transfer matrices, these being defined in terms of the bulk and boundary weights 
of the previous sections. Double-row transfer matrices of this type were introduced 



in and first used for interaction-round-a-face models in |33 . 

An additional feature of the lattice to be considered here is that alternate rows 
will be associated with spectral parameter values u and A— u. Also, the left boundary 
will be associated with the spectral parameter value A— u and the right boundary with 
value u. These values are used since they result in the double- row transfer matrices 
forming a commuting family. The most physically relevant point is the isotropic 
point u = A/2, at which all rows and both boundaries are associated with the same 
value of the spectral parameter. Combining this with the conformal point described 
in the previous section, the case of most interest here is thus u — £i = £2 = A/2. 

We begin by defining a set of paths consistent with boundary conditions (r l5 a\) 
and (r 2 , a 2 ) and lattice width N, 

^ = {(A, b ,..., b N , (3 2 ) I (60, • • • , b N ) G F^FZa, > °' 

{!,..., F- a j}. ^ 



We see that 



K iai \r 2a2 \ = (F ri G N F^) aia2 , (3.53) 
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which is invariant under interchange of a\ and a 2 or of r\ and r 2 , since F ri , F r ' 2 and 
G are symmetric and mutually commuting matrices. 

We now introduce a double-row transfer matrix D^ a A r2a ^u, £1, £ 2 ) with rows and 
columns labeled by the paths of Gr iai \r 2 a 2 anc ^ entries defined by 

2 (u, £1, 6)(/3i,6 ,...,6jv,/32),(5i,(io,...,djv,52) = 

d <5i 



riai| r2a2^ 



E 



5 



Co 



(c ,..., Cjv ) e 5"+' 

rAf 



A-u,fi x 



<n;=o« 



AT-1 

J"=0 



1) 



c i c i+i 



A — -u 



b N (3 2 



(3.54) 



5! dp i i i fio di 



djy-i d N d N g 9 



j n,oi \ 


\—u 

X 




A— u 

X , 


/ r 2 ,a 2 j 


\-u, £i/ 












u 




u 






X 




\ 





/3i b bo &i 



6jv-i b N b N 02 



We note that if £/£ 0l i r202 is empty, then 0^ ai | r2a2 (w, £i, £2) is zero-dimensional with 
its value taken as 0. 
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The partition function for the lattice model is now given by 



tr(^ ai|r2a2 Kei,e 2 )) J 





A— it 

\ 






A— u 

\ 


u 

N , 






N 










A— u 






A— u 


u 

\ , 






\ , 





2M (3.55) 



AT 



We can see from the first equality of (|3.55 ) that the task of evaluating the partition 
function is equivalent to that of evaluating the eigenvalues A^.^^w, £1, £2)*;, k = 
1, . . . , {F"G«F») aiaa , of Ju, £1, 6), with 



3.8 Transfer Matrix Properties 



Z ( A f iai [r 2 a 2 (w^l, 6)fc) 
fc 



(3.56) 



We shall now show that the double-row transfer matrices satisfy a variety of impor- 
tant properties. In particular, we shall find that applying certain transformations to 
the parameters of the model results only in similarity transformations of the transfer 
matrix. This implies that these parameter transformations are symmetries of the 
model, since the transfer matrix eigenvalues and partition function are invariant 
under any such similarity transformation. 



30 



3.8.1 Commutation 



It follows from the Yang-Baxter equation ( 3.12 ), inversion relation ( 3.13Q and bound- 
ary Yang-Baxter equation ( p. 33 ) that the double-row transfer matrices commute for 
any two values, u and v, of the spectral parameter, 



riai\r 2 a2^ 



o. 



(3.57) 



This can be proved diagrammatically, as done in Section 3.4 of |R3 



3.8.2 Transposition 



It follows straightforwardly from ( |3.10| ) and ( |3.38| ) that each double-row transfer 
matrix is similar to a symmetric matrix. More specifically, defining 



L r iai \r 2 a2 (f3 1 ,b ,...,b N ,p 2 ),(Si,do,...,dN,$2) 



1>. 



1/4 °(Pi,b ,--;b N ,l32),(S 1 ,d ,...,d N ,52) 
bo 



and 



(3.58) 



A-iai|r- 2 a 2 ( M '^l'^2) 



AN J~fN ( p p \ (AN \-l 

riai|r2a2 rioi|r2a2\ Sl> S2J \**Tiai\r202 ) ' 



(3.59) 



we have 



(A-iai|r 2 a 2 ( W '^l'^2) 



~ N 

D rxai\r 2 a2\ u iiliZ,V 



(3.60) 



This implies that if D na , r2a2 (u, £i, £2) is real, as for example occurs if ip is the 
Perron- Frobenius eigenvector and u, £1 and £ 2 are appropriately chosen, then the 
eigenvalues of ■Q^ 01 | r2a2 («, 61, £2) are & U rea l- 



3.8.3 Gauge Invar iance 



We see that the model is invariant under any gauge transformation ( 3.35Q of the 
boundary weights since this results only in a similarity transformation of the double- 
row transfer matrix, 



Dr iai \r2a 2 ( U -> £l> £2) ^ ^riai\r 2 a2 ^r 1 a 1 \r2a 2 ( U ^ £l> £2) (<S^ai| 



r2a2 / 



(3.61) 
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where 

CUV 

r iai \r 2 a 2 (/3 1 ,& ,...,6 iv ,/3 2 ),(<Si,do v ..,djV,<S2) ~ ^ g 2 ) 

S nai (&0W1 5 ,r2a2 (&Af)/3 2 <5 2 <W -,M ,(*,-, ^iv) • 

3.8.4 Simplification at Completely Anisotropic Points 



It follows from the anisotropy property ( |3.11| ), boundary anisotropy property ( |3.49| ) 



crossing symmetry ( |3.10|) and boundary crossing symmetry ( |3.39|) that at the com 



pletely anisotropic points, u = and u = A, the double-row transfer matrices are 
proportional to the identity, 

n jv (n c c \ r»Jv (\ c c\ ^ 2 s o(£i) s o(6) s ri(£l) s r- 2 (6) r / Q fiQ x 
Hiai[r 2 a 2 l U ) a;«J - M-iai|r 2 a 2 (A Sl> W - S (2fi) So( 2 6) 

3.8.5 Crossing Symmetry 

It follows from the Yang-Baxter equation ( 3.12 ), inversion relation (|3.13j ) and bound- 



ary crossing symmetry (|3.39|) that the double-row transfer matrices satisfy crossing 
symmetry, 

D r N iailr2a2 (\-u, 6, 6) = D r N iailr2a2 (u, Ct, 6) • (3.64) 
This, like (|3.57|) , can be proved diagrammatically, as done in Section 3.3 of |33| . 



3.8.6 Left-Right Symmetry 



It follows from reflection symmetry ( |3.9| ) and boundary reflection symmetry ( 3.38 ) 



that on interchanging the left and right boundary conditions and boundary fields, 
we have 



^r 2 a 2 |nai( W ' &j £l) ^f l0l |r 2 a 2 ~ ^riai|r 2 a 2 (-^nai|r 2 a 2 (^ ~ u i 6)) , (3.65) 

where i?^ ai | r202 is a square matrix with rows labeled by the paths of ^ 02 | riai , 
columns labeled by the paths of £/£ ai i r202 and entries given by 

R r iai \r 2 a 2 ( / 3 1) b 0) ...,6 N ,/3 2 ),(5 1 ,(io,...,d N ,5 2 ) = ^2,6jv,-,bo,/9i),(5i,do,-,diV,-52) ■ (3.66) 



32 



Combining ( |3.65| ) and the invertibility of R^ iai \ r2CL2 with Q3.60| ) and ( |3.64| ), we see that 
D r iai \r 2 a 2 ( u >Z,u&) and D r 2 a 2 \nai( u > 6, £1) are related by a similarity transformation, 
so that this complete interchanging of the left and right boundaries is a symmetry 
of the model. 



3.8.7 Symmetry under Interchange of r 1? £1 and r 2 , £2 

It can also be shown that under interchange of r\, £1 and r%, £2, we have 

A-20i|rio 2 (' U ' ^2, £l) Crlai|r 2 a 2 (6,6) = ^10^202(6 ' ^2) - D riai|r- 2 a 2 (' U ' &) , (3.67) 

where C^ ai | r2a2 (£i, £2) is a square matrix with rows labeled by the paths of ^ ai | n02 , 
columns labeled by the paths of Qr iai \r 2 a 2 an< ^ entries given by 



*-^riai|r 2 a 2 l 



(C C \ - ^ b Td Td N 
Ul ) ?2j(/3i,6 ,.-,6jv,/3 2 ),(<5i,rfo,-,rfiv,5 2 ) — ~ ~ x 



4>a? ^a 2 



f>0 c b N-l c N-l co d O c N-l d N-l 



E E •E E ••• E w 

(co,...,c W -l) 6 Q N ao=1 a *>-l =1 70=1 7jv_i=1 



r 2 ri 



/ do 7o c 
V ai /3i 6 



6-61 
+A 



X 



JV-1 

n 

i=o 



/ Cj 1 c i+ i 
\bj 1 6^1 







£ 2 +a1 




7j 7jfi 






\ c i 1 c i+i 





(3.68) 



in which we take 0^ = 02,0^ = /3 2 and 7^ = 5 2 - 

The weights which appear in (|3.68|) are fused bulk weights which are defined, for 
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each r, s £ {1, . . . , g— 1}, as fused r — 1 by s — 1 blocks of bulk weights, 
/ g? 7 c 



5 /3 



9i 



U s (a,d)s^ aMy ,„ )hg _ 24 ) 



Qr-2 



u+ 

(s-r)X 
% 




(s-2U 

s 








U— 

(r-2)A 




\ 


ei 


e, 


--2 



/s-2 



(3.69) 



^ s (6, c)p,{b,h,...,f^c) 



h 



U 7 (d, 6)a,(a,ei,...,e^2,6) 

These fused weights can be regarded as lattice generalizations of the 1 by r — 1 and 
r— 1 by 1 fused blocks considered in (2.20). They satisfy a generalized Yang-Baxter 



equation, as given in (3.39) of f4"4H , and it is by using this equation, and expressing 
the boundary weights in the double row transfer matrices in the form ( |3.32|) , that 
Q3.67D can be proved. In doing this, it is the (£i + 2A)-dependent 1 by r\ — 1 fused 
weights in ( |3.68| ) which, together with the generalized Yang-Baxter equation, allow 
the ^-dependent r% — 1 by 2 fused block of bulk weights within the left boundary 
weight to be propagated to the right. Similarly, it is the (£ 2 +A)-dependent 1 by r 2 — 1 
fused weights which allow the ^-dependent r 2 — 1 by 2 fused block within the right 
boundary weight to be propagated to the left, and it is the (£i — £ 2 + A)-dependent 
r 2 — 1 by r\ — 1 fused weight which allows the two blocks of bulk weights from the 
left and right boundary weights to be interchanged. 

We also note that the generalized Yang-Baxter equation allows the fused weights 
within C* iai \ r2a2 (£i, £2) to be rearranged. In particular, the (£1 — £ 2 + A)-dependent 



fused weight, which in ( |3.68| ) is on the far left of the lattice row, can be propagated 
to an arbitrary position further to the right, while reversing the order of the (£i+2A)- 
and (£ 2 +A)-dependent fused weights to the left of its final position. This essentially 
corresponds to the fact that, in proving ( |3.67|) , the order with which the blocks from 



34 



the left and right boundaries are interchanged is arbitrary. 

By using a generalized inversion relation, it can also be shown that C* iai \ r2a2 (£,i, £2) 
is invertible (except at isolated values of £1 or £ 2 which, by continuity, are unimpor- 
tant), its inverse being proportional to C r ^ ai | ria2 (£2, £1). It thus follows that transfer 
matrices with r±, £1 and r 2 , £2 interchanged are related by a similarity transforma- 
tion, so that this partial interchanging of the left and right boundaries is a symmetry 
of the model. 

Finally, we note that a simple generalization of the results of this section is that it 
is also possible, up to similarity transformation, to propagate the fused block of bulk 
weights within either boundary weight to an arbitrary position within the interior of 
the transfer matrix. The importance of this observation is that it is consistent with 
the viewpoint in conformal field theory of the boundary conditions corresponding 
to local operators. 

3.8.8 Properties Arising from £ — > ioo 

Various properties follow from the £ — > ioo (or equivalently £ — > —ioo) form (|3.50|) 
of the boundary weights. 

For example, applying this limit to both the left and right boundary fields we 
find that D^ ia ^ r2a2 (u,ioo,ioo) is proportional to a direct sum, over a{,a{ e Q, of 
D* ai n la3 i(u, £1, £2), with this term being repeated F^ ai ,F^ a2 , times in the sum. 

As a further example, we can attach an additional r\ — 1 by 2 fused block on 
the left and an additional r 2 — 1 by 2 fused block on the right of D^ l/a < r2 , a2 (u, £1, £ 2 ), 
where each of these blocks is of the same form as that in Q3.32j ) and we take £ — > ioo 
in each. We then find, from (|3.50|) , that the resulting matrix is proportional to a 
corresponding direct sum and, using methods similar to those of Section |3.8.7| , that 
the two attached blocks can be interchanged up to similarity transformation. This 
therefore gives the result 

0—'—' Kw\r 2 'aA^ £1, 6) « D r », ai , lr2 , a2 ,(u, £ 1; £ 2 ) , 

(3.70) 

where ~ indicates equality up to similarity transformation and the superscripts on 
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indicate the number of times that the corresponding terms appear in the direct 
sum. This relation implies that 



ai' ,ad&G 



^ Faiai' ^a 2 a 2 ' Zr 1 >"a 1 '\r 2 >a 2 '(. U i£,li&) 



(3.71) 



3.8.9 Decomposition into Single-Row Transfer Matrices at the 
Conformal Point 



It follows from ( |3.47| ) that if both boundaries are at their conformal point, that 
is if u = A — £i = £2 = then the double-row transfer matrix decomposes into a 
product of two single-row transfer matrices. We denote these, generally not square, 

V 

matrices as T* iai \ r2a2 (£), with rows labeled by the paths of Qr iai \r 2 a 2 anc ^ columns 

A 

labeled by the paths of G^+i ,ai|r 2 +i ,a 2 » 

and T ? iai \ r2a2 (£) , with rows labeled by the 
paths of ^+1,01^2+1,03 anc ^ columns labeled by the paths of Q, 
these matrices are given by 



n ai \r 2 a 2 - The entries of 



v 

rj-iN 



riai\r 2 a 2 



(0(ft 



b ,... ,6j V ,^ 2 ),(7i,c ,... ,cjv,72) 



N-l 



j=0 \°3 

0, f[G bjC] =0 

3=0 



n E r ^(b N ,c 



A 'ft 72) IT GbjCj 

3=0 



(3.72) 



7l r 



Co Cl 



cjv-i cjv 



I 
I 



ft b b a 



\ 



,72 



1 

nT2 , 12 I 
I 



fejv_i b N b N ft 
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and 



-^riai|r2a 2 (0(7liCo,->CN)72),(|8i,bo,->&JV>/82) 

bj b j+1 



JV-l 

^(^cow n W\ 

3=0 

N 

o, n^=o 

3=0 

/3i bp b p bi 



Cj Cj+i 



3=0 



I 

I 



71 



\ 



I 

T 2 , <J2 | 
I 



CQ CI 



cjv-i cat 



72 



(3.73) 



We now see, using ( |3.47| ), that the decomposition of the double-row transfer 
matrix at the conformal point is 



riai|r 2 a 2 l 



2 (£, A— £, £) — T riai]jJ , 2a2 (£) 2^ iai | r2a2 (0 • (3.74) 
We emphasise that the orientations of the boundary edge weights in the two single- 



row transfer matrices in ( |3.74|) are opposite with respect to a fixed direction along 
the boundary. Thus, even at the isotropic point £ = A/2, at which adjacent rows of 
the lattice become indistinguishable in the bulk, the alternating orientations of the 
boundary edge weights will still distinguish adjacent rows at the boundaries. 

We see from O that T^JC) = (A- +ljai|r2+lia J- 2 (T- ai|r2a2 (0) T (< ai | r2a2 ) 2 , 

where A^ ia ^ r2Cl2 is given by ( |3.58|) , and therefore that 

T 



"^riai|r 2 a 2 (^> > ^ ("^riaj r2« 2 ) "^riaijf^oaC^) ("^"ldil^"/ 



(0 (^riai|r 2 a 2 (0 J J ^r 1 a 1 \r 2 a 2 > (3.75) 



where T^,,^ 



/t\ _ /liV npjv fdUiV \-l 

: VSJ - f:L riai|r 2 a 2 ■ x riai|r 2 a 2 VS A nai|r 2 a.2 / 



This immediately implies that if T noi | r2a2 (£) is real, then -D£ 01 | raa2 (£, A-£,£) has 
all nonnegative eigenvalues. 
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3.8.10 Properties Arising from Graph Bicolorability 

We now consider some properties which arise if the graph Q is bicolorable. Although 
we have not assumed until now that Q is bicolorable, all of the specific A, D and E 
cases to be considered in the next section have this property. 

Bicolorability of Q means that a parity n a G { — 1,1} can be assigned to each 
node a G Q so that adjacent nodes always have opposite parity; that is, 

G ab = 1 =>• 7T a 7T b = -1 . (3.76) 

We now note that it follows from ( p . 1 8[ ) that each F r with r even/odd is an 
odd/even polynomial in G, which leads to a generalization of ( |3.76[ ) to a selection 
rule on the fused adjacency matrices, 

F r ab >0 n a 7f b = (-l) r+1 . (3.77) 

Proceeding to the effect of graph bicolorability on the lattice model, a square 
lattice can be naturally divided into two interpenetrating sublattices, with the near- 
est neighbors of each site on one sublattice being sites on the other. Thus, if Q is 
bicolorable then the condition that the spin states on neighboring lattice sites be 
adjacent nodes on Q implies that, in each spin assignment, the spin states on one 
sublattice all have parity 1, while those on the other all have parity —1. 

Using ( p.77| ), there is also a consistency condition between the left and right 
boundary conditions, 

= (-l) n+r * +N , (3.78) 

since otherwise </£ a]|raaa is empty and D^^u,^,^) and Z™ ]r2a2 (u, £i, &) are 
zero. If this condition is satisfied, then the sublattices are fixed by the boundary 
conditions; that is, the sublattice which contains all parity 1, or all parity —1, spin 
states is the same for all possible assignments. 

Finally, we note that ( |3.77| ) also implies that for each boundary edge (b, c) G £ m , 
the parities of b and c are respectively opposite to and the same as 7i a (— l) r , and 
that therefore any pair of nodes can appear in only one order in a particular set of 
boundary edges. 
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4 . Critical Unitary A-D-E Models 



In this section, we specialize to the critical unitary A-D-E models, for which Q is an 
A, D or E Dynkin diagram with Coxeter number g and ip is the Perron-Frobenius 
eigenvector of the adjacency matrix. In these cases, A = ir/g and the regime of 
interest here is < u < A. This class of models was first identified and studied in 
111 Oil. The A and D models are the critical limits of models introduced in IHTfl and 



respectively, but the E models do not have off-critical counterparts. We also 
note that if tp is instead taken as a submaximal adjacency matrix eigenvector, then 
the critical nonunitary A-D-E models are obtained. 



Fusion of the A models was introduced in 41, 42, 43 1, while fusion of the D and E 



models was first studied in ||44|| . We shall also consider intertwiner relations among 
these models, these having been studied in detail in [^3|, [4^, |51]]. In particular, we 
shall find that certain symmetries in the fusion and intertwiner properties of these 
models lead to various additional properties of the boundary conditions, transfer 
matrices and partition functions. 

We shall also study the explicit forms of the boundary weights and boundary 
edge weights for these models. For the A models, various methods have previously 
been used to obtain sets of diagonal boundary weights in 
containing non-diagonal weights in M, 31, 



32, 33, M and sets 



and all of the A boundary 
weights found here represent certain cases of these previously-known weights. For 
the D and E models, sets of diagonal boundary weights were found by direct solution 
of the boundary Yang-Baxter equation in |53] |, but most of the sets obtained here 
contain non-diagonal weights and were not previously known. 

Finally, we shall consider the connection between the lattice model boundary 
conditions at the conformal and isotropic points and the conformal boundary con- 
ditions of the corresponding unitary minimal theories. 
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4.1 Graphs and Adjacency Matrices 

The A, D and E Dynkin diagrams with Coxeter number g are explicitly given by 



9-2 9-1 ' 



+1 



■2.-2 

2 z 2 1 




5 = 2,3,... , 



9 = 6,8,... 



(4.1) 



(4.2) 



3 + 
2 



1 2 3 4 5 ' 
7 



E 7 



1 2 3 4 5 6 ' 



9 = 12, 



9 = 18, 



(4.3) 



(4.4) 



and 



Ea 



9 = 30. 



(4.5) 



3 4 5 6 7 ' 

We note that when referring to the nodes |~ and | + of we shall use the 

convention that if the superscript is not specified then the choice is immaterial. 



Furthermore, if a numerical value for |~ or | is required in any equation, it is to 
be taken as |. 

The eigenvalues of the adjacency matrices of these graphs are 2 cos(kj n/g), j = 



,\Q\, where kj are the Coxeter exponents, as given explicitly by 



A 



9-1 



k i = J 



kj — 



'2j 



1 ?' = 1 2 



. 2 ' •> 



i?6 

Ea 



(h,...,h) = (1,4,5,7,8,11) 

(h,...,k 7 ) = (1,5,7,9,11,13,17) 

(ki, ...,k s ) = (1,7,11,13,17,19,23,29). 



(4.6) 
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It can also be shown that for an integer k, and any A-D-E graph with Coxeter 
number g, 

3 an adjacency matrix eigenvector with all non- 

_ . , . . •<=>- k is coprime to o . (4.7) 

zero entries and associated eigenvalue zcosykn/g) 

We note, as already mentioned in Section |3.4.1| , that ([4.7D implies that, for any 
A-D-E case, the maximum fusion level and Coxeter number are equal, and that 
they can thus both be denoted by g. 

Throughout the rest of Section 4, we shall restrict our attention to the unitary 
A-D-E models, which are associated with the Perron- Frobenius eigenvector and 
thus with the Coxeter exponent k% = 1. Therefore, from this point on, we shall take 

A = n/g (4.8) 

and if) to be the Perron-Frobenius eigenvector, which can be explicitly given by 



2=1,.. 



.4 



9-1 



2-1 
' 2 x 



Da. 



+i 



a 
a 
a 



1,... 

1-2 

l-l 



1-3 



(4.9) 



Ei, I 



6,7,8 



S a , a = 1, 
,l/(2sinA), a = f 

sin((/-l)A)/sin(2A) 
2cos((/-2)A) , 

i sin((/-3)A) /sin(2A) , a = l 

where, from ( |2.4| ), S a = sin (a A) / sinA. 

It is also important here to consider a particular involution a i— > a of the nodes 
of each graph, this being given by the graph's Z 2 symmetry transformation for the 
A, -Dodd and E e cases and by the identity for the D even , E 7 and E 8 cases. More 
explicitly, we have 

a E Ag_i 



g-a 

2 ' 

6 — a 
a . 



e Di 



2 ' f+1 J 

a e E 6 \{6} 
otherwise . 



f + 1 odd 



(4.10) 
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We see that the eigenvector entries fl4.9|) are invariant under this involution, 

i>a = i>a. (4.11) 

We also observe that each A, D and E graph is bicolorable and that we may set 
the parities as 

ir a = (-l) a . (4.12) 

We now consider the A, D and E fused adjacency matrices. A more comprehen- 
sive treatment of these can be found in Appendix B of 0. 
For Ag-i, we have explicitly 



ab 



(4.13) 



1; a+b+r odd, \a — b\<r — 1 and r + 1 < a + b < 2g— r — 1 
; otherwise . 

We note that F£ b for A g -i, with r ^ g, is actually symmetric in all three indices. 
For Do -i-i, we have explicitly 



'2; a, 6^1, a + 6+r odd, |a — 
and a + > lf-rl + 1 



b\<l 



1; a,6^f, a + 6+r odd, |a - 6| < f -|f-r| -1 
and |a+6-|| < ||-r|-l 



ab 



1; a^f, 6 



2 ' 



a+^+r odd and a > || — rl+1 



(4.14) 



|, 6+| + r odd and 6>||-r| + l 



I*, & = f ± and r = 1 mod 4 



2 ' 



| T and r = 3 mod 4 



1 ; a 
1 ; a 
1 ; a 

; otherwise . 

The fused adjacency matrices for E§ are given explicitly in Section [4.4.6| . We 
shall not give the E>? and E$ fused adjacency explicitly, since they can be obtained 
straightforwardly using a computer from the recursive definition ( |3.18| ), but we note 
that for E 7 each entry is in {0, . . . , 4} and that for E 8 each entry is in {0, . . . , 6}. 

We now list some properties which apply to all of the A, D and E fused adjacency 
matrices. These properties can be proved by decomposing these matrices in terms 
of their eigenvalues and eigenvectors as given in M. 
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For r = g we have 

F 9 = 0. (4.15) 

Meanwhile, for r£ {1, . . . , g— 1}, F r form the basis of a commutative matrix algebra, 
which is a representation of the Verlinde fusion algebra. These g — 1 matrices are 
therefore often referred to as Verlinde matrices and denoted V r . 

These matrices can also be used to define related intertwiner matrices I a , for 
each a e Q, with rows labeled by 1, ... , g—1, columns labeled by the nodes of Q and 
entries given by 

I a r b = K b - (4-16) 

Each I a then intertwines the fused adjacency matrices of Q and those of 

i a F(gy = F(A g ^yr. (4.17) 

For Q = Ag-i, this property simply amounts to the commutation of the A g _\ fused 
adjacency matrices, but if Q is a D or E graph it forms the basis of various rela- 
tionships between Q and the A graph with the same Coxeter number. 

Finally, a property of the A, D and E fused adjacency matrices of particular 
relevance here is that, for r e {1, . . . , g — 1}, 

Fy = Fr ab = F a V (4.18) 

An important special case of this is r = 1, for which 

= S ai . (4.19) 

4.2 Bulk Weights, Fusion Matrices and Fusion Vectors 

We now consider the bulk weights, fusion matrices and fusion vectors of the critical 
unitary A, D and E models. 
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We see from ( |3.6| ) and ( |4.9|) that the A g _i bulk weights are 
ail a 



W\ 

w\ 



a a=Fl 

a ail 
a=|=l a 

a ail 
ail a 



si{-u) 

(Sa-! S a+1 )^ S (U) 



(4.20) 



s a (±u) 

S a 



It can also be shown, using (|2.15| ), (|3.8| ) and ( |4.20|) together with certain results 
on the fusion of bulk weights from |[43|| , that the A s _i fusion matrices are given 
explicitly by 



P r ( a , 0)(a,ci,...,c r _ 2 ,6),(a,di,...,d r _ 2 ,fe) 



r-j-a— 6— 1 
2 



r-a+b-l 
2 



r-1 



n m n ^ m 

m=2 / V m=2 / 



n s m n 



r-1 



_ a + i)-r+l 



^m) 1 ^, 



n ^ 



(4.21) 



. m=2 



where in the fourth product we set Co = d = a and c r _i = d r -i = 6, and where 

(4.22) 



1, a = or 1 mod 4 
— 1, a = 2 or 3 mod 4 . 



In fact, the only properties of the sign factors required here are e a G { — 1, 1} and 
£a-i e<n-i — ~1> so an y of the three other assignments which satisfy these could be 
used instead. 

We can see from ( [4.21 ) that, in keeping with ( |jj.2(J| ) and ( 4.13| ), each nonzero 
A g _i fusion matrix has rank 1 and that the corresponding fusion vectors, which are 
thus uniquely defined up to sign, are given, up to this sign, by 



U T (a, 0)l,(a,ci,...,c r _2,fe) 

r+a— b— 1 

n s m )( 



m=2 



i — a + b-l 

n s r , 

m=2 



a+h+r-l 

2 



r-1 



n 5 « 

. a+b-r+l 



n s - 



m=2 



1/2 / r-1 



n 



(4.23) 



nift / ' 
, m=0 " Jc "t / 



where in the last product we set Co = a and c r _i = b. 

Proceeding to -Df+i, it is possible to write expressions, similarly explicit to those 
for Ag-i, for the bulk weights and entries of the fusion matrices and fusion vectors, 
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for a certain natural choice of these vectors. Since each of these expressions involves 
many different cases, analogous to those of ( |4.14j ), we do give them here. However, 



we note that, as with A g _i, each Dz +1 bulk weight and fusion vector entry can be 
expressed as a single product of terms. 

We also note that each Ds +1 bulk weight and fusion matrix entry can be ex- 
pressed as a linear combination of A g -\ bulk weights and fusion matrix entries, with 
the coefficients being products of entries of so-called intertwiner cells. These inter- 
twiner cells, whose exact properties are outlined in [49], |50, 51fl , serve a similar 



role at the level of the bulk weights to that served at the level of the adjacency ma- 
trices by the intertwiner matrices ( |4.16|) . While the expressions provided by these 
intertwiner cells may not be particularly compact when applied to specific cases, 
they are still useful for deriving certain general properties using known properties 
of the intertwiner cells and of the A g ^\ bulk weights and fusion matrices. 

With regard to the E graphs, the number of different cases is particularly large 
so that the numerical values of fusion matrix and fusion vector entries are probably 
best evaluated using a computer. While this may result in a somewhat unnatural 
choice of the fusion vectors, this is largely immaterial since the lattice properties of 
interest are independent of this choice. We also note that, exactly as with -D|+i, 
each E bulk weight and fusion matrix entry can be expressed as a linear combination 
of Ag^i bulk weights or fusion matrix entries using intertwiner cells and that these 
expressions can be used to obtain general properties of the E bulk weights, fusion 
matrices and fusion vectors. 

We now consider some properties which apply to all of the A, D and E graphs. 
We first note, using (|4.11|) , that the bulk weights are invariant under the involution 





u) . (4.24) 



We also see similarly that for the fusion matrices, 

-P r («, &)(a,ci,..., Cr _2,fc),(Mi,...,d r -2,b) - P r (^^)(a,c 1 ,...,c r -2,b),(a,d 1 ,...,d r - 2 ,b) ■ ( 4 -25) 

It follows from this that U r (a, &) Qi (a,ci,...,c r _ 2 ,fe) an d U r (a, b) a ,(a,ci,...,cr-a,i) correspond to 
two orthonormal decompositions of both P r (a,b) and P r (a,b), and therefore that 
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these fusion vectors are related by a transformation ( |3.25| ), 



U r (a,b) a ^ £l ^- Cr _ 2! - b) = ^2 K r (a,b) aa , U r (a,b) Q ,^ Cl ^ Cr _ 2:b) . (4.26) 

a'=l 



Finally, we note that the fusion matrices P 9 1 (a,a), which from (|4.19| ) all have 
rank 1, can be used to generate the fusion matrices of lower fusion level according 
to 



E 



P 9 (a, a.)(a,ci,...,c r _2,b,ei,...,e g _ r _2,a),(a,di,...,d r _2,6,ei,...,e 9 _ r _2,a) 



Sri> a 



P r (a, b) 



(4.27) 



(a,ci , . . . ,c r _ 2 ,t>) , (a,di , . . .,ti r _2 ,£>) 



This can be proved by first obtaining the result for the A graphs using ( 4.21|) and 
then proceeding to the D and E graphs using intertwiner cells. 



4.3 Additional Boundary Condition and Transfer Matrix 
Properties 

We now show that in addition to satisfying all of the properties outlined in Sec- 
tion [3.8| , including those of Section |3.8.10| arising from bicolorability with parity 
fl4.12 ), the critical unitary A-D-E models possess further important symmetries 
associated with the involutions ( [4.101) and the intertwiner relations ( |4.16|) . 



4.3.1 Symmetry under a,\ i— > a x and a 2 



a 2 



We first consider the relationship between the (r, a) and (r, a) boundary conditions. 
It follows straightforwardly from ([4.18|) and (|3.41|) that the sets of boundary edges 
for these boundary conditions are related by 



{{b,c)\(b,c)eS ra } 



(4.28) 



Proceeding to the corresponding boundary edge weights and boundary weights, 
we find, using ( ^4.26|) , that 



E™(b,c) 



/3'=1 7 '=1 



(4.29) 
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and 



B ra \ c r 



d 6 



b 13 



)^ba ^da / 
= E E S ra (b)^, S ra (d) S5 , B ra (c 
/3'=1 <5'=1 \ 



d 5' 

b (3' 



(4.30) 



where the orthonormal transformation matrices in ( (4.26 ) and those in ( 4.29|) and 
( ^301) are related by (|3~37|) . 

It now follows from fl424|) and flUD that D^-Ju, 6, 6) and D^Ju, 6) 
are related by a similarity transformation, 



D N (11 P P \ I / S 1 ^ 

"^riai|r2a2\ ' >1> riai\r20,2 riai|r2Q2 



rj-iV ^iv t-jJV (up P \ 

riai\r2d2 riai|r2a2 "^riai|r2ti2\ ' >2J 



(4.31) 



where S^ iai{r2a2 is given by (|3.62Q , using the same S ra (e) as in ( |4.30| ), and H^ ai{r2a2 
is a square matrix with rows labeled by the paths of £/^ 5l i r2S2 , columns labeled by 
the paths of ^ 0l | r802 and entries given by 



if 



riai|r 2 a 2 (fofio,...^ ,f3 2 ),(Sl,do,-,d N ,h) 



5, 



(/3i,6 0l — ,6jv,/9i),(5i,do,— 4n,$2) 



(4.32) 



4.3.2 Equivalence of (r, a) and (<?— r— l,a) Boundary Conditions at the 
Conformal Point 

We now show that, at the conformal point and for r G {1, . . . , g— 2}, the (r, a) and 
r— l,a) boundary conditions are equivalent. This equivalence takes the form of 

the boundary edge weights for each of these boundary conditions being the same, 

except for a gauge transformation and a reversal of the order of the nodes in each 

boundary pair, neither of which affects any properties of interest. We shall denote 

such an equivalence of boundary conditions by <->. 

In terms of the set of boundary edges ( |3.41| ), it follows from ( |4.18| ) that, for 

r G {1, . . . , g— 2}, the sets £ ra and S 9 ~ r ~ 1,a contain the same pairs of nodes but with 

opposite ordering; that is, 



S g-r~xra = {( Cj6 )|( 6>c ) e s ra }. 



We also see, from ( ^.15| ), that 



£ 



g-l,a 



(4.33) 



(4.34) 
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so that at the conformal point there are no (g— 1, a) boundary conditions. 

For the boundary edge weights, the equivalence of the (r,a) and (g — r — l,a) 
boundary conditions is 



E 9 



— r— 1,5 



F ba Ka 1 

[c,bU = E E S ra (b)^, S r+1 > a (c), Y E™(b,c)^,, 
p'=\ 7 '=i 



(4.35) 



where S ra (e) are orthonormal matrices satisfying ( 3. 36 ) which are defined by 

S ra (e) aa , = (5 r VaM) 1/2 x 

E E U9 ~ («) a )l,(o 1 6 g ^-2,...,6l,e,ci,..., C ^ 2 ,a) X (4.36) 

(e,fei,...,fe 9 _ T ._ 2 ,a) egf,r r (e,ci,...,(V-2,o) 6 Q% a 



We note that S ra (e) also satisfy 

S ra {e) T = S 9 - r ~ a {e). 



(4.37) 



These relations, ( 4.35 ) and ( (4.37] ), and the ortho normality ( p. 36 ) can all be proved 
using the general properties of the fusion matrices and fusion vectors, (|3.15| ), ( |3.16| ), 
( PT2TD , ( ^22j) and those which follow from ( PT7D , together with ( |Q7D . 

We now note that, as labels, (r, a) and r— l,a) are always distinct since for 
an A odd , D or E graph g is even so that r and g—r—1 are different, while for an A even 
graph g is odd so that a and a = g — a are different. Thus, due to the equivalence 



( 4.35 ), there are at most (g — 2) |£/|/2 distinct boundary conditions. Furthermore, 
by examining the sets of boundary edges (|3.41 ) for each A, D and E case, we find 
that there are no further equivalences between these sets, either direct or through 
reversing the order of nodes in each edge. We therefore conclude that 

number of boundary conditions at the conformal point = (g — 2) \G\/2 . (4.38) 



Due to the consistency condition 



, the implementation of a given left and 



right boundary condition at their conformal points on a lattice of fixed width can 
be achieved using only two of the four possibilities which arise from the two versions 
of each boundary condition. If (n, ai) and (r 2 , a?) is one of these possibilities, then 
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the other is (g— T\ — 1,Oi) and (g — r2 — 1,^2) and the transfer matrices for the two 
are related by 



g—ri—l,ai\g—r2—l,a2\ 



(4.39) 



-i,, 2 (£,A-£,0 = 

- \ A v JV 1 

^ri+l,oi|r2+l,02 T^ iai \ r2a2 (\— T^ ia ^ r2a2 (\— £) (<Sv. 1+1)ai | r . 2+1)02 ) , 

where S r ^ +ljai | r2+lj(l2 is given by ( p . 6 2|) and (|4.36| ). Comparing this with ( |3.74|) we see 
that the ordering in the products of single row transfer matrices is different in each, 
which does not affect the nonzero eigenvalues (although unimportant differences 
in the number of zero eigenvalues will arise due to the different dimensions of the 
oppositely-ordered products). We thus see that the partition functions are related 
by 

Z^-iM^-id^ W, = z^rJl-t W) • ( 4 - 4 °) 

Finally, we find, using ( [2.30|) in (|3.28|) , that the boundary edge weight relation 
(|4.35[) implies the boundary weight relation 



d 5 

b (3 



«,d = E E s ra (b) Pf3 , s ra {d) ssl b™ (c 



/ f3'=l 8'=1 

with S ra (e) again given by ( |4.36|) . This then implies that 



d 5' 

b (3' 



(4.41) 



' riai\r2d2 



' g— ri,ai| r2a2 



1 riai\g—r2,a,2 



(4.42) 



This relation taken at u — A — £1 = £2 = A — £ is consistent with ( |4.40|) through the 
equivalence, which follows from ( |3.48|) , of the (r,a) boundary condition at u = — £ 
and the (r — 1, a) boundary condition at u — £. 



4.3.3 Intertwiner Symmetry 

We now consider the relationship between the double-row transfer matrices and 
partition functions of the model based on Q, with Coxeter number g, and those of 
the model based on A g _i. In particular, we shall find that any critical unitary A— 
D-E partition function can be expressed as a sum of certain A partition functions. 
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It can be shown using intertwiner cells that, for each ri,r 2 , s' G {1, . . . ,g— 1} 
and ai, a 2 G Q, we have 



9-1 



F - («, a, 6) « («, a, 6) 



8=1 

9-1 



(4.43) 



8=1 



where we are using pa and the superscripts on in the same ways as in (|3.70|) , and 
where the fused adjacency matrices on both sides and the transfer matrices on the 
left sides refer to Q, while the transfer matrices on the right sides refer, as indicated, 

tO Ag-\. 

Since the two forms of this relation can be proved similarly, and are related 
through the symmetries of Sections 3.8.6 and 3.8.7] , we shall consider, from now on, 
only the first form. 

We first discuss the details of the proof. The similarity transformation in the 

first line of ( [4.43 ) can be implemented by an invertible matrix J s r ^j r2a2 which pre- 

multiplies the left side and postmultiplies the right side. The rows of J^^nas anc ^ 

9-1 nA 

the rows and columns of F ^ D^ s ,^ (u, £1, £ 2 ) are labeled by the paths of 

s=l 

{(to, ...,t N ,t, a 2 ) I (1, t , ••• , t N , 1) G ( Vi)? lS 'M> F L 2 > 0, a 2 G {1, . . . , F* ia j} , 

while the columns of ^'101^202 anc ^ tne rows anc ^ c °l umns of F " ia D* a < r2a2 (u, £1, £ 2 ) 
are labeled by the paths of 



{(Cti, 6 , • • (Pl,bo, . . . ,b N ,/3 2 )eG? lb \r 2a2 , F a' 1 b>®> 0!iG{1, 

The entries of these matrices are given by 



nai|r 2 a 2 (t ,...,t N ,t,a 2 ),(a 1 ,b,l3 1 ,b ,...,b N ,/3 2 ) 



a l b 


a l b N 


( b fa 6 \ 


N-l 


( bj 


bj+i \ 




^ &jv /?2 «2 \ 


E ■ 


. . Q riai 


ai 70 


UQ ai 




7i+i 


Q r2ai 


7iv "2 


70=1 


7iv=l 


{s' 1 t J 


j=0 


\tj 


tj+i / _ 




V** 1 t J 



(AAA) 
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Laeg J (ai,6,/9 1 ,6 ,...,6 N ,/3 2 ),(a! / 1 ,d,(5 1 ,(io,- 1 ^,52) (4.45) 



J OL\OL\ 



and similarly for Faia2 i^^2sH u 5 £i) £2)- In (|4.44j) , Q ai are intertwiner cells as- 

s=i 



sociated with the intertwiner matrix I ai of ( 4.16| ), and Q Tiai and Q r2ai are fused 
blocks of such cells, of widths r\ — 1 and r 2 — 1 respectively. Thus, J^^^ can be 
viewed as a row of intertwiner cells in which the spins on the lower left and upper 
right corners are fixed to s' and a 2 respectively, the lower row of spins between s' 
and a 2 on A g _i label the rows of the matrix and the upper row of spins between s' 
and a 2 on Q label the columns of the matrix. We also see that all of the matrices 
here are square with dimension (I ai F ri G N F r2 ) s > a2 , it being possible using ( [4.171 ) to 
propagate I ai to the right of this expression while replacing the adjacency matrices 
of Q with those of A g -\. The intertwiner cells are assumed to satisfy an intertwiner 
relation, as given in (4.6a) of as well as two inversion relations, as given in 
(4.6b) and (4.6c) of |[49|| . It follows immediately from the inversion relations that 
Jr 1 a 1 \r 2 a 2 * s invertible, with its inverse being given, up to gauge transformations on 
the intertwiner cells, by its transpose. Meanwhile, the equation corresponding to the 
first line of (|4.43|) can be obtained by using one of the inversion relations to insert 
a pair of cells between and F " ia i^ a | r2a2 (w, £i, 6), using the intertwiner 

relation and the form ( |3.32| ) of the boundary weights to propagate these cells around 
a single loop, and then using an inversion relation again to remove the inserted cells, 

thus giving F " ia2 D^ s ,^{u, £i, £ 2 ) Jr&ifaca anci completing the proof. We note 

s=l 

that this last process is probably best understood diagrammatically, with the lower 
part of the loop involving the conversion in the lower row of the transfer matrix 
from Q weights to A g _i weights and the propagation of the row of intertwiner cells 
to a position between the rows of the transfer matrix, and the upper part of the 
loop involving the conversion of weights in the upper row of the transfer matrix and 
the propagation of the row of cells to its final position. 



We note that (|4.43|) is still nontrivial for Q = A g _\ and that it then corresponds 
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to certain cases of fl3.70p . In fact, the A g _i intertwiner cells Q a can be obtained by 
taking a u — > ioo limit on fused a — 1 by 1 blocks of A g -\ bulk weights. 

We also note that the intertwiner cells Q a for Q a D or E graph have only been 
found explicitly, in |23|, f|9], |50], pTJ, for a = 1. However, for various reasons, the 



existence of these cells for other values of a seems guaranteed. 

Finally, we observe that a particularly important case of fl4.43|) is s' — 1, for 
which 

D^\r 2 aM^& « © F -'« ^afcWi.GO , (4-46) 
which in turn implies that 

K^u, a, 6) = E < a2 («, 6, &) • (4-47) 

8=1 

We thus see that the task of evaluating the partition functions of all of the A-D-E 
models with boundary conditions (r\,ai) and fa, 0,2) has been reduced to that of 
evaluating the partition functions of just the A models with boundary conditions 
(ri,l) and (r 2 ,s). 



4.4 Boundary Weights 

We now consider the explicit forms of the boundary weights and boundary edge 
weights for the critical unitary A-D-E models. 

4.4.1 A Graphs 



For Ag-i, we find, using Q4.23Q in (|3.42j ) and then applying a simple gauge transfor- 



mation to remove a factor e c which arises, that the boundary edge weights are given 
explicitly by 

ET(c±l,c) u = (^c+a)/2 S^ y2 )» 

V ; (S c± iS c )^ v ; 

We see that each of these weights is positive. We also see that, for these weights, 
the relations (|42g ) and flP5p are 



E r *- a {g-b,g-c)n = E ra (b,c) u , E^ r -^- a {c,b) xl = E ra (b,c) u . (4.49) 
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B 



Substituting ( |4.48| ) into ( |3.46| ) we now find that the A g _i boundary weights are 

v c±1 1 , 

S( r ^c+a)/2 S( c ± aT r)f2 Sq(£,+u) S r (£ — u) + S 



(r±c+a)/2 S( CT a±r)f2 



B' 



c±l 1 



S r (S c S c±1 )^s (2£) 

-c+a)/2 S( r +c-a)/2 5'(c+a-r)/2 S ( c+a+r )/2) 1/2 Sq(2u) 



(4.50) 



(S c ^S c+1 )^S^s (2C) 
All of the boundary weights which were found as solutions of the boundary Yang- 
Baxter equation for the A g -\ models and their off-critical extensions in |4|, [31], |32|, |33|, 



34] , p5| can be related to those of ([4.50|) by using appropriate values for the various 
parameters involved. In particular, we note that more general boundary weights 
which depend on two boundary field parameters are known, as for example given 
in (3.14-3.15) of 0], and that these reduce to the weights (|4.50|) when one of these 
parameters is set to zoo. 

Some important special cases here are the (r, 1) and (g — r — l,g — l) boundary 
conditions for which 

S rl = {(r,r+l)}, = {( r+ i, r )}, 

E rl (r,r + 1) = Es- r - 1 >e-i(r+l,r) = (S r S r+1 ) 1/4 . 
The corresponding boundary weights are all diagonal and given by 



(4.51) 



B rL r±l 



r 1 
r 1 



B g-r,g-ll r±1 



r 1 
r 1 



u. 



S 1 ^ s (£±u) s r (Z,Tu) 
Sf *„(20 



(4.52) 



these weights matching those found in p3| . We see that for r G {2, ...,g — 2} 
these cases provide an example in which boundary weights which are nonzero away 
from the conformal point vanish at the conformal point. More specifically, away 
from the conformal point, these cases represent semi-fixed boundary conditions in 
which the state of every alternate boundary spin is fixed to be r and that of each 
spin between these can be r — 1 or r + 1, while at the conformal point they rep- 
resent completely fixed boundary conditions in which only a single boundary spin 
configuration . . . r, r + 1, r, r + 1 . . . contributes to the partition function. 
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Other important cases here are the (l,a) <-» (g — 2,g — a) boundary conditions, 
with a G {2, . . . , g — 2}, for which 

£ lQ = {( a ,a-l),{a,a+l)} , £9-2,9-a = {( _i ; a ) ; ( a+ i ; a )} ; 

(4.53) 

E la (a,a±l) = E°- 2 *- a (a±l,a) = (S a±1 /S a ) 1/4 . 

As already seen for the general case in (|3.31|) , these are semi- fixed boundary condi- 
tions in which the state of every alternate boundary spin is fixed to be a. 

We now observe that the only cases in which every edge of A g -i appears, in one 
order, in the set of boundary edges are 



g odd: 



g even: 



= {(1,2), (3, 2),..., (0-2, 0-3), (0-2, 0-1)} 
£^i,4i = {(2,1), (2,3), ...,(0-3,0-2), (0-1, 0-2)} 

E^A = {(2,1), (2, 3),..., (0-2, 0-3), (0-2, 0-1)} 
£M = {(l,2),(3,2),...,(0-3,0-2),(0-l,0-2)}. 



(4.54) 



These cases therefore represent boundary conditions in which each configuration of 
boundary spins consistent with fixed even-spin and odd-spin sublattices contributes 
a nonzero weight to the partition function at the conformal point. If these weights 
are all equal, which in fact only occurs for A3, we refer to the boundary condition 
as free, while if the weights are not all equal we refer to it as quasi-free. 

Finally, we note that all other boundary conditions not of the fixed, semi-fixed, 
quasi-free or free type can be regarded as intermediate between these types. 

4.4.2 Example: A 3 

We now consider, as an example, the case A 3 , this being the simplest nontrivial A 
model. For any spin assignment in the A 3 model, the spin states on one sublattice are 
all 2, while each spin state on the other is either 1 or 3. The former sublattice, being 
frozen in a single configuration, can therefore be discarded and the model viewed as 
a two-state model on the other sublattice. It can be shown that the bulk weights 
of this model are those of the critical Ising model, with the horizontal and vertical 
coupling constants suitably parameterized in terms of the spectral parameter. It is 
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therefore natural to relabel the nodes of A 3 as a frozen state and Ising states + 
and — ; that is, 

A 3 = • . . . (4-55) 



The A 3 fused adjacency matrices are, from ( 4.13|) 








l A = I I . (4.56) 




Using these matrices to determine each set of boundary edges and then (|4.48|) to 
give the corresponding weights, we find that these weights are, up to normalization, 
as given in Table 1. 





E(-,0)n = 1 


£(0,+)n = 1 





B(0,+)n =£(0-)n = 1 


£(+,0)u =£(-,0)11 = 1 


+ 


£(+,0)n = 1 


E(0,-)u = 1 


a / 

/ r 


1 


2 



Table 1: A 3 boundary edge weights 

We see that the three A3 boundary conditions at the conformal point are: 

• (1,+) <-> (2,-) <-> + fixed 



;i,o) 



;2,o) 



— fixed 
free. 



We note that the last of these boundary conditions is invariant under the model's 
Z2 symmetry, while the first two are not. 



4.4.3 Example: A 4 

We now consider, as another example, the case A4. It is known that this model can 
related to the tricritical hard square and tricritical Ising models. 
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The fused adjacency matrices are 



^ 1 ^ 
10 
10 



/ 1 ^ 
10 10 
10 1 



F 3 



\0 1 / 


\ 1 


J \ 




/o 1 ^ 




/oooo\ 


F 4 = 


10 


F 5 = 







10 









V 1 J 




v 0/ 



/OOIO^ 
10 1 
10 10 
10 0/ 



(4.57) 



Using these and ( 4.48|) , we find that the A 4 boundary edge weights are, up to 
normalization, as given in Table 2. 



4 


£(4,3)n = 1 


£(3,2)n = 1 


£(2,1)11 = 1 


3 


£(3,2)n = (v^+l) 1 '* 
£(3,4)n = (v^-l) 1 -* 


£(2,l)n =£(4,3)n 
= (V5+1)V8 

£(2,3)n = (5 v / 5-ll) 1 '* 


£(l,2)n = (v^-l) 1 / 8 
£(3,2)n = (v^+l) 1 / 8 


2 


£(2,l)n = (VS-l) 1 ^ 
£(2,3)n = (v^+l) 1 ^ 


£(l,2)n =£(3,4)n 

= (v^+l) 1 / 8 
£(3,2)n = (Sv^-ll) 1 ^ 


£(2,3)n = (v^+l) 1 / 8 
£(4,3)n = (VE-l) 1 / 8 


1 


£(l,2)n = 1 


£(2,3)n = 1 


£(3,4)n = 1 


a / 

/ r 


1 


2 


3 



Table 2: A 4 boundary edge weights 
We see that the six A4 boundary conditions at the conformal point are: 



• (1,1) - 


+ (3,4) <- 


...1,2,1,2 


. . fixed 


. (1,4) «- 


- (3,1) <- 


...3,4,3,4 


. . fixed 


. (2,1) <- 


+ (2,4) «-» 


...2,3,2,3 


. . fixed 


. (1,2) <- 


- (3,3) <-> 


...2,1/3,2, 


1/3 . . . semi- fixed 


. (1,3) <- 


- (3,2) <- 


...3,2/4,3, 


2/4 .. . semi- fixed 


» (2,2) <- 


+ (2,3) <-> 


quasi-free. 
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4.4.4 D Graphs 



For Ds +1 , we find, by substituting the fusion vector entries described in Section [O 
into ( |3. 42 ) and then applying certain simple gauge transformations, that the bound- 
ary edge weights can be taken explicitly as 



E ra (b, c) 



11 



f 2 -l/4 
2 -l/4 

2 1/4 



6,c=f 



6,c=f 



EX a (b,c) u - 



< 2-1 



c)n; a^f, r>§ 



J .4 

9 



b,C=f 



£4 2 (Mil 



2 
2 



£ ra (6, c) 



12 



fO; Ml 

±2 -i/4 6 



r < z 

2 ' ' — 2 



±2 -i/4 



11 



2± r > 9 
2 ' ' — 2 



(4.5£ 



c) 21 



fo; c ^f 

±2 -i/4 ^-«fa_ 1 r 



^4 V 2 _ x ' 2^ n ' C 



±2 -i/4 £ ;m ( |_ 1) | )ii 



4± r <2_l 

2 ' ' — 2 x 

r > 3- 

2 ' ' — 2 



£ ra (&, c) 22 



11 



r<f-l 



,i?l a (&,c) 



11 



> i 



2 • 



where by X\ bc= g we mean that X is only to be included if b = § or c 



|, and where 



E™(b, c)n are the boundary edge weights as given by ( f4.48| ). The fact that 
the -D|+i boundary edge weights can be expressed in terms of A g _ 1 boundary edge 
weights follows from the intertwiner relations between the corresponding models. 

We see that for any boundary condition, and in this gauge, there is at most one 
negative boundary edge weight. We also see that, for these weights, the relations 
( gggD and flP5|) are 

E r ~ a (b,c)^ = a Pl E ra (b,c) Pl , E^-^icb)^ = E ra (b,c)^, (4.59) 



where 



a 9l 



- 1, D odd with (3 ^ 7 
1, otherwise. 



(4.60) 
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In fact, for -D ev en the first relation of ( |4.59| ) is trivial since the involution a i— ► a is the 



identity, but we note that the relation still holds if the involution is instead taken 
to be the graph's Z 2 symmetry transformation and if a 12 — °"2i are replaced by — 1. 

Boundary weights for -D2+1 can be obtained by substituting the edge weights 
Q4.58|) into (|3.46|) . For some of the boundary conditions, all of the boundary weights 
are diagonal and these weights can be related to those previously found in . How- 
ever, non-diagonal boundary weights for the D& +1 models, apart from one case of D4 
considered in |7j], were not previously known. We also note that since some simple, 
but (r, a)-dependent, gauge transformations have been included in the boundary 
edge weights of ( |4.58| ), some corresponding gauge factors need to be included in 
equations which relate boundary weights at different values of (r, a) , in particular 



Finally, we note that, as with A g -i, certain of the -D|+i boundary conditions 
can be identified as being of fixed, semi-fixed, free or quasi-free type. 

4.4.5 Example: _D 4 

We now consider, as an example, the case -D4. For any spin assignment in the 
D4 model, the spin states on one sublattice are all 2, while each spin state on the 
other is 1, 3 or 4, so that the model can be regarded as a three-state model on the 
latter sublattice. The bulk weights of this model can be shown to be those of the 
critical three-state Potts model. These bulk weights are also invariant under any 5 3 
permutation of the Potts spin states. We shall use the more natural labeling of the 
nodes of -D4, 1 1— > A, 2 1— > 0, 3 + 1— > B and 3~ 1— > C; that is, 



( ggp and (1441) . 




(4.61) 
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The -D4 fused adjacency matrices are, from ( 4.14|) 



F 5 



/ 1 0\ 
10 
10 

Vo 1 J 



/o 1 0\ 
10 11 
10 

Vo 1 0/ 



(4.62) 



F 3 



/o 1 1 ^ 




( \ 


2 


F 5 = 





10 1 







V 1 1 J 




v / 



where the rows and columns are ordered A, 0, B, C. Using these and ( |4.58| ), we 
find that the boundary edge weights are, up to normalization and a simple gauge 
transformation on the (2, 0) and (3, 0) boundary conditions which makes their £3 
symmetry properties more apparent, as given in Table 3. 



c 


E(C,0) u = 1 


#(0,^)11 = 
E(0,B) u = 1 


E(A,0)u = 
E(B,0)u = 1 


E(0,C)n = 1 


B 


E(B,0)u = 1 


E(0,A) n = 
E(0,C)u = 1 


E(A,0)u = 
E(Cfi)u = 1 


E(0,B)u = 1 





E(0,A) u = 
E(0,B)n = 
E(0,C)n = 1 


E(A,0)u = 1 
E(A,0) 12 = 
E(B,0)n = -1/2 
E(B,Q) 12 = V3/2 
E(C,0) n = -1/2 
E{C,0) 12 = -V3/2 


E(0,A)n = l 
E(0,A) 2 i = 
E(0,B)u = -1/2 
£(0,5)31 = V3/2 
E(0,C) u = -1/2 

£(0,(7)21 = -V3/2 


E(A,0)u = 
E(B,0)n = 
£(C,0)n = 1 


A 


E(A,0)u = 1 


E(0,B) n = 
£(0,0)1! = 1 


E(B,0)u = 
E(C,0)n = 1 


£(0,A)u = 1 


a / 

/ r 


1 


2 


3 


4 



Table 3: -D 4 boundary edge weights 
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We see that the eight D4 boundary conditions at the conformal point are: 



(M) 


<— > 


(4, A) 


<-> A fixed 


(1,B) 


<— > 


(4,5) 


<-> B fixed 


(hC) 


<— > 


(4,C) 


<-> C fixed 


(2, A) 


<— > 


(3, A) 


<-> i? and C mixed with equal weight 


(2,B) 


<— > 


(3,B) 


<-> A and C mixed with equal weight 


(2,C) 


<— > 


(3,C) 


<-> A and S mixed with equal weight 


(1,0) 


<— > 


(4,0) 


<->• free 


(2,0) 


<— > 


(3,0) 


<-> quasi-free in which same-spin pairs 



The nature of the last boundary condition is best seen by considering the (2, 0) 
boundary weights at the conformal point, these being, up to normalization, 



This is therefore a boundary condition on nearest-neighbor pairs of Potts spins, in 
which like and unlike neighbors are associated with weights 1 and —1/2 respectively. 

We see that the last two -D4 boundary conditions are S3 symmetric, while the first 
six are not. In fact, the (2, 0) boundary weights (|4.63|) represent the only possibility, 
other than reproducing the (1,0) «-» (4,0) weights, which is S3 symmetric and 
consistent with a decomposition ( |3.47|) in which 7 is summed over two values. 

4.4.6 E Graphs 

For the E graphs, there are a large number of boundary conditions at the conformal 
point (specifically, 30 for E$, 56 for E-j and 112 for Eg) and many of these are in 
turn associated with a large number of boundary edge weights. Therefore, since it is 
straightforward and more practical to obtain the numerical values for these weights 
using a computer, we do not list them here. However, as an example, we do give 
the sets of Eq boundary edges. 



and different-spin pairs have weight —1/2 





(4.63) 
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The E 6 fused adjacency matrices are 



F 1 = 



F 4 = 



F 7 = 



/ 1 n 


n 

u 


n 

u 


n 

u 






/ n 


i 


n 
u 


n 

u 


n 
u 








n 

u 


1 

J. 


n 
u 


n n \ 




1 
















1 





1 


















1 





1 


1 




n n 


1 

1 


n 
u 


n 
u 


n 


F 2 = 


n 


1 

1 


n 
u 


1 

1 


n 
u 


1 

1 


F 3 = 




1 


n 
u 


9 


n 
u 






n n 


u 


1 


n 

u 


n 




n 


n 
u 


1 

1 


n 
u 


1 

1 


n 

u 







1 


n 
u 


1 

1 


n 1 













1 
















1 


















1 










\0 











ly 




\o 





1 








o) 




VO 


1 





1 


J 




/ n n 


n 

u 


1 

1 


n 
u 


i \ 

L 




/ n 


n 
u 


1 

1 


n 
u 


1 

1 


n \ 




/o 


1 

1 


n 
u 


1 

1 


n n \ 







2 





1 










1 





2 





1 






1 





2 





1 




n 9 

u z 


n 

u 



z 


n 

u 


i 

_L 


F 5 = 


i 


n 

u 


Q 
o 


n 

u 


1 

1 


n 

u 


F 6 = 







9 
z 


n 

u 


9 
z 


n 9 

u z 




1 n 



z 


n 

u 


n 

u 


n 




n 


9 


n 

u 


1 

1 


n 

u 


1 

1 






1 


n 

u 


9 
z 


n 

u 


i n 




1 











1 




1 





1 


















1 





1 







\i o 


1 





1 


Oy 




\o 


1 





1 





1/ 




Vo 





2 





o o y 




/l 


1 








ON 




/ n 


1 

1 


n 

u 


n 


n 

u 






f° 





1 





0\ 




2 





1 





1 




1 





2 


















1 





1 


1 




1 


3 





1 





F 8 = 


n 


9 


n 

V_7 


9 


n 


1 

J. 


F 9 = 




1 





2 





1 




1 





2 





1 


n 
u 


n 




Z 




i 


n 







1 





1 


1 







1 





1 







o 


o 


o 


1 


o 


1 












1 










Vo 1 





1 





ly 




U 





1 





1 


Oy 1 




Vo 


1 





1 


J 




/ 1 \ 




/o 











1 


o\ 






/000000\ 





1 





1 
















1 



























1 





1 





1 


F n = 








1 











F 12 = 



















1 


1 
















1 

































1 
















1 




































V 1 J 




Vo 














u 






Vooooooy 



(4.64) 



Using these, we find that the E 6 boundary edges are as given in Table 4. In this 
table we also give, for each (6, c) G £ ra , the values of F fe r a and F^ 1 as successive 
superscripts. 
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6 


(6,3)" 


(3,2)" 
(3,4)" 


(2,1)" 
(2,3)" 
(4,3)" 
(4,5)" 


(1,2)" 
(3,2)" 
(3,4)" 
(3,6)" 
(5,4)" 


(2,3) 12 
(4,3) 12 
(6,3) 12 


(3,2) 21 
(3,4) 21 
(3,6) 21 


(2,1)" 
(2,3)" 
(4,3)" 
(4,5)" 
(6,3)" 


(1,2)" 
(3,2)" 
(3,4)" 
(5,4)" 


(2,3)" 
(4,3)" 


(3,6)" 


5 


(5,4)" 


(4,3)" 


(3,2)" 
(3,6)" 


(2,1)" 
(2,3)" 
(6,3)" 


(1,2)" 
(3,2)" 
(3,4)" 


(2,3)" 
(4,3)" 
(4,5)" 


(3,4)" 
(3,6)" 
(5,4)" 


(4,3)" 
(6,3)" 


(3,2)" 


(2,1)" 


4 


(4,3)" 
(4,5)" 


(3,2)" 
(3,4)" 
(3,6)" 
(5,4)" 


(2,1)" 
(2,3) 12 
(4,3) 12 
(6,3) 12 


(1,2) 12 
(3,2) 22 
(3,4) 21 
(3,6) 21 


(2,1) 21 
(2,3) 22 
(4,3) 12 
(4,5)" 
(6,3) 12 


(1,2)" 
(3,2) 21 
(3,4) 22 
(3,6) 21 
(5,4) 12 


(2,3) 12 
(4,3) 22 
(4,5) 21 
(6,3) 12 


(3,2) 21 
(3,4) 21 
(3,6) 21 
(5,4)" 


(2,1)" 
(2,3)" 
(4,3)" 
(6,3)" 


(1,2)" 
(3,2)" 


3 


(3,2)" 
(3,4)" 
(3,6)" 


(2,1)" 
(2,3) 12 
(4,3) 12 
(4,5)" 
(6,3) 12 


(1,2) 12 
(3,2) 22 
(3,4) 22 
(3,6) 21 
(5,4) 12 


(2,1) 21 
(2,3) 23 
(4,3) 23 
(4,5) 21 
(6,3) 13 


(1,2) 12 
(3,2) 32 
(3,4) 32 
(3,6) 32 
(5,4) 12 


(2,1) 21 
(2,3) 23 
(4,3) 23 
(4,5) 21 
(6,3) 23 


(1,2) 12 
(3,2) 32 
(3,4) 32 
(3,6) 31 
(5,4) 12 


(2,1) 21 
(2,3) 22 
(4,3) 22 
(4,5) 21 
(6,3) 12 


(1,2)" 
(3,2) 21 
(3,4) 21 
(3,6) 21 
(5,4)" 


(2,3)" 
(4,3)" 
(6,3)" 


2 


(2,1)" 
(2,3)" 


(1,2)" 
(3,2)" 
(3,4)" 
(3,6)" 


(2,3) 12 
(4,3) 12 
(4,5)" 
(6,3) 12 


(3,2) 21 
(3,4) 22 
(3,6) 21 
(5,4) 12 


(2,1)" 
(2,3) 12 
(4,3) 22 
(4,5) 21 
(6,3) 12 


(1,2) 12 
(3,2) 22 
(3,4) 21 
(3,6) 21 
(5,4)" 


(2,1) 21 
(2,3) 22 
(4,3) 12 
(6,3) 12 


(1,2)" 
(3,2) 21 
(3,4) 21 
(3,6) 21 


(2,3)" 
(4,3)" 
(4,5)" 
(6,3)" 


(3,4)" 
(5,4)" 


1 


(1,2)" 


(2,3)" 


(3,4)" 
(3,6)" 


(4,3)" 
(4,5)" 
(6,3)" 


(3,2)" 
(3,4)" 
(5,4)" 


(2,1)" 
(2,3)" 
(4,3)" 


(1,2)" 
(3,2)" 
(3,6)" 


(2,3)" 
(6,3)" 


(3,4)" 


(4,5)" 


a/ 

/ r 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 



Table 4: E 6 boundary edges 



From this table, the properties Q4.2S ) and ( 4.33 ) are immediately apparent, and we 
can also gain some understanding of the nature of each of the 30 boundary conditions 
at the conformal point. 
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4.5 Realization of Conformal Boundary Conditions 

We now consider the relationship between the integrable boundary conditions of the 
critical unitary A-D-E lattice models and the conformal boundary conditions of 
the critical series of s£(2) unitary minimal conformal field theories. 

Each conformal field theory of this type on a torus or cylinder is associated with 
two graphs, the A graph A g _2 and an A, D or E graph Q with Coxeter number g, this 
theory being denoted Ai(A g _ 2 , Q). As shown in [jTBl [T7|, [TJ|, the lattice model based 
on Q, with if) the Perron- Frobenius eigenvector and < u < A, can be associated 
with the field theory A4(A g _2, G)- 

In 0, it was found that the complete set of conformal boundary conditions of 
M.[A g -2i G) on a cylinder is labeled by the set of pairs (r, a), with r £ {1, . . . , g — 2}, 
a e G and (r, o) and (g — r — 1, a) being considered as equivalent. We immediately 
see that this classification is identical to that of the boundary conditions of the 
corresponding lattice model at the conformal point. 

It was also shown in || that the partition function of Ai(A g _ 2 , G) on a cylinder 
with conformal boundary conditions (r 1; ai) and (r 2 ,a2) is given by 

Zr^naM = E E f(VaL nG) s aia2 X(r,.)(?) , (4-65) 

r=l s=l 

where g is the modular parameter, F(y4 9 _ 2 ) r and denote the fused adjacency 

matrices of A g _2 and G, and X(r,s) is the character of the irreducible representation 
with highest weight 

(rg -s(g-l)) 2 
of the Virasoro algebra with central charge 



A (r,s) = 7T7I i\ (4.66) 



c = 1-^-rr- (4-67) 
0(0-1) 

The equivalence of the (r, a) and (gr— r — l,o) conformal boundary conditions is 
apparent by using ( |4.18| ) and the relation A( g _ r _ l 9 „ s ) = A( r s ) to observe that the 
partition function (|4.65|) is unchanged by applying this transformation to either of 
the boundary condition labels. 
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We now assert that, in the continuum scaling limit, the (r, a) boundary condition 
in the lattice model at the conformal and isotropic point provides a realization 
of the (r, a) conformal boundary condition of the corresponding field theory. In 
particular, we expect that, as iV — > oo, the eigenvalues of D* ia < r2a2 (\/2, A/2, A/2) can 
be arranged in towers, with each tower labeled by a pair (r, s) and the multiplicity of 
tower (r,s) given by F(A g -2) r r r F{Q) s aia2 - We further expect that the j'th largest 
eigenvalue in this tower has the form 

AjjjU = exp[-2^-2/^ + Mc/24-A (r> . ) -* (r ,.y)/^ + o(l/JV)] l (4.68) 

where T is the bulk free energy per lattice face, which depends only on g, f ri \ r2 is 
the boundary free energy per lattice row, which depends only on g, r\ and r 2 , c and 
A( rjS ) are as given in (|4.66|) and (|4.67|) , and k(r >s )j are nonnegative integers given 



through the expansion of the Virasoro characters by 



oo 



X(r,.)(9) = q~ C/24+A ^ £ 9* (r - W , hr,s)j < hr,s)j+t ■ (4-69) 

3=1 

With the eigenvalues appearing in this tower structure, it follows using (|3.56|) 
that, as N — > oo, the lattice model and conformal partition functions are related by 

^WA/2, A/2, A/2) ~ 

exp(-2MW-2M/ ri | r2 ) Z riai , r2a2 (g), g = exp(-2vrM/Ar) . 

We note that we also expect related conformal behavior in a lattice model which is 
at the conformal point and has < u < A, but which is no longer at the isotropic 
point u = A/2. 

The expectation that the lattice model boundary conditions correspond in this 
way to the conformal boundary conditions is supported by the results of numerical 
studies we have performed, the matching of the identifications made here of the 
nature of the lattice realizations of certain conformal boundary conditions with 
those made in other studies, the consistency of all of the symmetry properties of the 
lattice model partition function with those of the conformal partition function, and 
analytic confirmation in several cases. 
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In our numerical studies, we evaluated the eigenvalues of D^ r2S (X/2, A/2, A/2) 
for certain A graphs, selected values of ri, r 2 and s, and several successive values of 



N. We then extrapolated these results to large N and verified consistency with ( 4.68 
for these cases. This numerical data, used with ( |4.46| ), also implied consistency with 
( f4.68| ) for all of the related A, D and E cases. 



Regarding the identification of the nature of the lattice realizations of particular 
conformal boundary conditions, this was done for all A 3 cases and all -D4 cases except 



(2,0) in 0, for all (l,a) and (r, 1) cases of A g ^ and D 2+1 in JZTj, for all A 4 



cases in ||25|| , and for the (2, 0) case of D± in | 29| . In all of these studies, the lattice 



model boundary conditions were shown to have exactly the same basic features as 
those found here. 

Proceeding to the consistency of symmetry properties, it follows straightfor- 



wardly from (|4.65| ) and the properties of the fused adjacency matrices that the 
Ai(Ag_ 2 ,Q) partition function satisfies 

^riai\r2a.2 \Q) ^g— ri— 1,01)5— T2— 1,52 (?) 
E F{SYa\ ai > F{G)Za 2 > Z r , aArM {q) = 

ai',a 2 'e6 (4-71) 

E F(gy a \ a , z ri , aiW ( q ) 

J2 F(g)i ia Z ria{r2a2 (q) = E F (GYaia2 Z ri s>\r2s{q) ■ 
aGG s=l 

We immediately see that these equalities are consistent with the lattice relations 
(HD, (ET67D , (|43lD , (|440| ), ( gjlD and (|443|) respectively. We also note that certain 
cases of the last equality of fl4.71|) and its lattice version ( [4.43D , mostly with s' = 
r i — r 2 — 1; were considered numerically and analytically in [[[], ^(| ^T], |22], 23, |24j| , 



while in free combinations of (1, a) lattice boundary conditions in A g _\ were 
studied analytically leading to ( J4.47| ) with r\ — r 2 = 1 and a sum on a x and a 2 . 
However, in all of these studies, the orientation of the lattice differed from that used 
here by a rotation of 45 degrees. 

Finally, the partition function relation ( }4.70| ) has been proved analytically using 
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techniques based on the Yang-Baxter and boundary Yang-Baxter equations for a 
lattice with the same orientation as that used here for all A3 cases in [^H] and for 
all A 4 cases which lead to a single character on the right side of ( |4.65| ) in JSTJ. 



5 . Discussion 

We have obtained various results on integrable boundary conditions for general 
graph-based and, in particular, critical unitary A-D-E lattice models. More specif- 
ically, we have systematically constructed boundary weights, derived general sym- 
metry properties and, for the A-D-E cases, studied the relationship with conformal 
boundary conditions. 

The general formalism presented here, or certain natural extensions of it, can be 
applied to various other integrable lattice models which are associated with rational 
conformal field theories of interest. We expect that the corresponding integrable 
boundary conditions provide realizations of the conformal boundary conditions of 
these theories, although we acknowledge that for many of these models only the 
bulk weights are currently known and that explicitly obtaining boundary weights 
may involve certain technical challenges. Nevertheless, in conclusion, we list these 
other cases and indicate their connections with those studied here: 

• If, using ( |4.7D , we choose an eigenvector of an A-D-E adjacency matrix corre- 
sponding to a Coxeter exponent k with 1 < k < g — 1 and k coprime to g, then a 
nonunitary A-D-E model is obtained which corresponds to the nonunitary mini- 
mal theory JH(A g -k-i, G)- This enables the consideration of all of the A4(Ah-i, G) 
minimal theories with h < g. The A4(Ah-i,G) theories with h > g and G a D or 
E graph can not be obtained in this way, but all of the Ai(Ah_i, A g _i) theories are 
accessible since in this case h and g are interchangeable. 

• By taking, in Section 3, G as an A^ l \ D^ l > or E^ 1 ' Dynkin diagram, and ip as the 
Perron- Frobenius eigenvector of its adjacency matrix, lattice models are obtained 
which correspond to certain conformal field theories with central charge c = 1. As 
noted in Section |3.4.1| , for these models s is given by the second case of Q2.2|) and 



there are infinitely many fusion levels. 
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• By replacing the relations of (|2.1 ) with those of the Hecke algebra, certain lattice 
models and conformal field theories associated with s£(n), for n > 2, can be ob- 
tained. In this case, although the lattice models become significantly different, we 
expect that most of the results of at least Section 2 remain unchanged 

• By using the dilute Temperley-Lieb algebra instead of the Temperley-Lieb al- 
gebra, lattice models which correspond to the so-called tricritical series of unitary 
minimal theories, Ai(A g ,Q), can be obtained. However, we note that the dilute 
Temperley-Lieb algebra contains considerably more generators than the Temperley- 
Lieb algebra, so that the formalism would be more complicated from the outset. 

• It is also apparent that lattice models whose bulk weights are given by fused 
square blocks of A-D-E bulk weights could be considered by applying some rela- 
tively straightforward extensions to the formalism. The field theories associated with 
these models include certain superconformal theories. The cases of the A models 



involving only diagonal boundary weights were studied in ||33|| , where it was shown 
that the double-row transfer matrices of the standard and fused model together 
satisfy a hierarchy of functional equations. We expect that equations of a similar 
form can be derived, using the boundary inversion relation ( |3.4(J ) and other local 
properties, for the remaining A-D-E cases, such equations being important in the 
analytic determination of transfer matrix eigenvalues. 

• Finally, we mention the off-critical A and D models, which can be associated with 
perturbed conformal field theories. The bulk weights in these cases can no longer be 
expressed in terms of the Temperley-Lieb algebra, but a fusion procedure still exists 
and boundary weights constructed from fused blocks of bulk weights attached to 
diagonal boundary weights still satisfy the boundary Yang-Baxter equation. Some 



integrable boundary weights are known for these listed in ||34j| , and it is 

expected that each of the A and D boundary conditions found here corresponds to 
a critical limit of an off-critical integrable boundary condition. 
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